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SECTION  I 
INTRODUCTION 


High  strain  rate  characterization  of  several  materials  has 
been  described  in  an  earlier  report  [1],  The  results  in  this 
report  include:  (1)  the  use  of  induction  coil  technique  to  perform 
high  temperature  split  Hopkinson  bar  tests,  (?)  the  use  of  manganin 
stress  gauge  technique  to  measure  the  stress  history  in  impact 
experiments,  (3)  an  experimental  technique  involving  a  flyer  plate 
impacting  a  long  rod  to  measure  initial  yield  stress,  ( J0  the 
measurement  of  shear  strength  using  transverse  and  longitudinal 
manganin  gauges,  (5)  an  improved  experimental  setup  for  expanding 
ring  experiments,  (6)  the  high  strain  rate  modeling  of  several 
metals  using  the  Bodner-Par tom  model  under  both  room  and  high 
temperatures,  (7)  the  precursor  decay  analysis  based  on  a  new 
offset  >  1  e 1 d  criterion,  and  (8)  a  new  dynamic  failure  model  Dasad 
on  an  improved  void  growth  model  approach. 

The  various  high  strain  rate  experimental  techniques  and 
results  for  several  metals  are  reported  in  Section  II.  Analytical 
modeling  of  the  high  strain  rate  data  is  reported  in  Section  111 
and  IV.  Summary  of  the  results  and  several  recommendations  are 
described  in  Section  V. 
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SECTION  II 
EXPERIMENTS 


For  high  strain  rate  material,  characterization,  we  employed 

the  split  Hopkinson  bar  (SHB)  technique  to  deform  specimen  under 

uniaxial  stress  in  both  tension  and  compression  at  strain  rates  of 
2  3 

<0  -10  /sec,  the  expanding  ring  technique  to  measure  the  average 
stress  obtained  at  strain  rates  of  above  5000/sec  and  at  large 
strains,  the  plate  impact  tests  to  measure  the  initial  yield  stress 
in  uniaxial  strain  at  very  high  strain  rates,  and  the  flyer  plate 
on  long  rod  impact  experiments  to  obtain  initial  yield  stress  at 
strain  rates  of  1 00 -500 /sec .  Manganin  gauge  techniques  to  measure 
stresses  in  the  plate  and  long  rod  impact  tests  are  discussed.  The 
use  of  the  induction  coil  heating  method  in  SHB  experiments  is  also 
reported.  These  techniques  have  been  described  in  detail  in 
References  1-4,  For  completeness,  a  brief  description  of  each 
technique  and  data  are  presented  in  this  section. 

2.1  SPLIT  HOPKINSON  BAR 

2.1.1  Basic  Configuration 

The  SHB  apparatus  consists  of  a  striker  bar  that  is 
made  to  impact  end-to-end  onto  the  pressure  bar.  The  impact 
produces  a  stress  wave  in  the  pressure  bar  whose  duration  is  twice 
the  acoustic  transit  time  of  the  striker  bar.  The  specimen  is 
placed  at  the  other  end  of  the  pressure  bar  as  shown  in  Figure  1. 
The  specimen  is  also  connected  to  the  transmission  bar.  The 
striker,  pressure,  and  transmission  bars  are  all  the  same  material 
and  diameter  (12.7  mm  diameter  Inconel  71 8).  'l’he  stress  wave  that 
travels  down  the  pressure  bar  is  partially  transmitted  and  par¬ 
tially  reflected  by  the  specimen.  The  reflected  and  transmitted 
waves  iri  the  bars  are  detected  by  strain  gages.  Analysis  of  the 
strain-gage  signals  yields  the  load  and  displacement  history  of  the 
specimen.  Detailed  discussions  on  the  strain-gage  analysis  can  be 
found  in  Refs,  5  and  6. 
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Figure  1.  Schematic  Sketch  of  the  SHB  Setup. 


For  tensile  measurements,  a  collar  is  placed  around 
the  threaded-end  specimen,  and  the  specimen  is  screwed  into  the 
pressure  and  transmission  bars.  The  initial  compression  wave  is 
transmitted  through  the  collar,  and  reflects  as  a  tensile  wave. 

When  the  tensile  wave  returns  to  the  specimen,  the  collar  falls 
away  and  only  the  specimen  sustains  the  tensile  load.  The  data 
from  the  Hopkinson  bar  strain  gages  are  used  to  calculate  engineer¬ 
ing  stress  and  strain  in  the  sample,  which  is  converted  to  true 
stress  and  strain  until  the  onset  of  localized  deformation 
(necking)  by  using  the  conventional  relationships. 

2.1.2  High  Speed  Photographic  Technique 

A  photographic  system  was  used  to  record  images  of 
the  specimen  profile  throughout  the  test.  The  system  is  described 
in  detail  in  Ref.  7.  Light-emi tting-diodes  (LED)  back  illuminate 
the  specimen  which  is  viewed  by  a  35  mm  rotating  drum  camera.  One 
photograph  of  the  necking  specimen  could  be  obtained  in  every  20 
microseconds  with  a  framing  rate  of  50  kHz.  The  Bridgman 
parameters  a  and  R  were  determined  from  the  various  photographs  at 
different  time  intervals  of  the  specimen  (see  Reference  8), 

Bridgman  equations  were  then  used  to  calculate  the  s tress- strai n 
relations  beyond  the  onset  of  necking. 

2.1.3  Induction  Coil  Heati ng  Techni que 

To  obtain  SHB  data  under  high  temperature,  an  induc¬ 
tion  coil  heating  technique  was  developed.  During  the  past  20 
years  many  workers  designed  various  heating  systems  in  order  to 
perform  hi gh- temperature  tests  in  the  SHB  system.  Heating  of  the 
specimen  was  achieved  by  placing  various  kinds  of  heaters  near  it. 
There  are  two  major  difficulties  with  those  techniques:  (1)  the 
time  needed  to  neat  the  specimen  is  relatively  long;  and  (2)  the 
heating  is  not  localized  and  a  large  volume  (including  bar  tips)  is 
heated.  Thus,  a  temperature  gradient  is  established  along  the 
loading  bars  to  either  side  of  the  specimen.  Such  gradients 
produce  variations  in  the  elastic  modulus  of  the  bars  and  hence 
variations  in  the  elastic  velocity.  As  a  result,  elastic  pulses 
propagating  along  the  bar  suffer  from  partial  reflections  and  from 
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change  of  amplitude.  Several  techniques  to  reduce  these  tempera¬ 
ture  gradients  and  methods  to  account  for  them  in  the  data 
reduction  process  are  described  in  Ref.  9.  The  need  to  heat  the 
specimen  at  high  rates  also  emerges  when  one  wants  to  avoid 
mi cros tructur al  changes  (like  grain  growth  and  recrystallization) 
which  can  occur  at  elevated  temperatures.  Rosenberg  et  al.  L 4  J 
developed  a  simple  heating  technique  which  is  both  fast  and 
localized . 

The  specimen  to  be  heated  is  surrounded  by  an  induc¬ 
tion  coil  which  carries  a  high-frequency  current  as  shown  in  Figure 
2.  The  specimens  were  heated  by  a  radio  frequency  generator 
manufactured  by  L ind berg/Cycle  Dyne  (Model  No.  A-50)  .  This  gener¬ 
ator  uses  the  240-V  line  power,  which  is  stepped  up  via  a 
transformer  to  about  5000  Vac.  This  high  voltage  is  rectified  and 
the  resulting  dc  is  supplied  to  an  RF  power  oscillator  operating  at 
approximately  450  kHz.  The  output  of  this  oscillator  is  about  5  k’W 
and  this  is  coupled  to  the  sample  by  using  0.25-in.  copper  tubing 
coils,  which  surround  the  sample.  The  output  power  is  usually  set- 
to  maximum  in  order  to  make  the  heating  process  as  short  as 
possi ble . 

Since  the  magnetic  field  varies  inversely  with  the 
square  of  the  distance  between  the  specimen  and  the  coil,  one 
wishes  to  work  with  the  smallest  possible  coll  if  the  heating  time 
is  a  limiting  factor.  The  coils  used  in  the  present  study  were 
larger  than  the  Hopkinson  bar  diameter  (12.7  mm)  by  only  about  3  rum. 
In  order  to  avoid  heading  the  collar,  we  replaced  the  metallic 
collar  by  a  ceramic  one.  Alumina  AD-995  (manufactured  by  Coors) 
was  used  because  it  has  an  acoustic  impedance  which  is  very  close 
to  that  of  the  Inconel  718  bars.  This  matching  of  impedance 
results  in  a  smooth  passage  of  the  3tress  wave  from  the  first  bar 
to  the  second  through  the  ceramic  collar. 

The  temperature  of  the  specimen  just  prior  to  loading 
was  measured  by  use  of  a  thermocouple  inserted  through  a  hole 
drilled  in  the  ceramic  collar.  The  thermocouple  touched  the 
specimen  so  that  the  temperature  was  monitored  continuously  during 
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the  heating  process.  When  the  temperature  reached  the  desired 
value,  the  striker  bar  was  fired.  In  order  not  to  heat  the  bars 
too  much,  we  heated  the  specimens  for  about  one  to  two  minutes  in 
all  the  experiments  described  here.  We  could  heat  steel  specimens 
to  700°C  in  about  one  minute  with  the  double  wound  arrangement.  We 
also  tried  a  single-coil  configuration,  but  this  was  not  as  effi¬ 
cient  as  the  double  coil  probably  because  of  worse  matching  between 
specimen  and  field.  In  some  of  the  experiments,  the  ends  of  the 
bars  were  heated  but  their  temperature  was  much  less  than  that  of 
the  specimens  and  the  heated  portion  of  the  bar  extended  to  only 
about  10  mm.  To  determine  the  effect  of  such  heating  on  the  wave 
propagation  characteristics,  we  conducted  two  experiments  in  which 
the  tips  of  the  bars  touched  each  other  and  were  heated  to  200°C 
and  il00°C .  The  transmitted  and  reflected  pulses  were  compared  with 
those  obtained  at  room  temperature  and  the  difference  in  the 
measured  strain  pulses  were  less  than  one  percent.  Since  the  SUB 
is  made  of  Inconel,  its  elastic  modulus  is  relatively  insensitive 
to  temperature  -  decreasing  by  only  20  percent  when  heated  to 
7Q0°C  . 

The  thermocouple  to  measure  the  temperature  in  the 
specimen  was  an  ungrounded  probe  (manufactured  by  Omega,  Probe  No. 
SCASS-03236) .  This  probe  is  especially  suited  because  it  can  touch 
the  specimen  without  causing  triggering  of  the  scope.  Other 
probes,  which  were  grounded,  caused  trigger  problems  because  the 
specimen  and  the  bar  are  part  of  the  triggering  circuit.  The 
maximum  reliable  temperature  with  this  probe  is  927°C,  according  to 
the  manufacturer.  Also,  due  to  the  very  small  dimensions  of  the 
probe  tip  (diameter  of  0.9  mm)  temperature  equilibrium  with  the 
specimen  was  fast.  In  order  to  estimate  the  uniformity  of  tne 
temperature  along  the  specimen  during  heating,  a  few  test  were  run 
in  which  the  temperature  was  checked  at  a  few  points  along  the 
specimen.  At  an  average  temperature  of  500°C  the  maximum  dif¬ 
ference  between  the  points  was  less  than  10°C.  A  value  of  ±  2 
percent  correction  may  be  needed  for  the  temperature  determination. 
We  conducted  high  temperature  SHB  tests  in  tension  on  several 
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metals.  Figure  3  summarizes  the  results  for  W-2  tungsten,  52125 
steel,  C1008  steel.  Nickel  200,  and  6061-T6  aluminum. 

2 . 1 . 4  Results 

Using  the  SHB  experimental  techniques,  compression, 
tension,  and  necking  tests  under  room  temperatures  and  tension 
tests  at  elevated  temperatures  were  conducted  on  several  metals. 
Table  1  described  the  various  tests  that  were  conducted  under  the 
present  study.  Quasi-static  tests  were  also  conducted  for  evaluat¬ 
ing  the  rate  dependency  in  these  materials  by  comparison  with  SHB 
data . 


a .  Armco  Iron 

Several  tension  and  compression  tests  were  con¬ 
ducted  at  room  temperature.  This  material  exhibited  a  fair  amount 
of  rate  dependency  as  can  be  seen  from  the  true-stress  vs.  log  e 
plot  in  Figure  H .  The  true  stress  vs.  true  strain  curves  at 
various  strain  rates  are  also  shown  in  Figure  5.  The  flow  stress 
under  SHB  strain  rate  levels  varied  between  6  and  7  kbars.  The 
strain  hardening  was  almost  negligible.  Due  to  the  strain  rate 
sensitivity  in  this  material,  the  SHB  flow  stresses  are  more  than 
double  those  of  the  quasi- static  tests. 

b.  Nickel  200 

Two  quasi-static  tests  at  strain  rates  of  0.001 
and  1.0s  1  and  several  SHB  tests  under  tension  were  conducted  at 
room  temperature.  The  results  are  shown  in  Figure  6.  The  stress- 
strain  curves  at  these  strain  levels  were  very  similar  as  can  be 
seen  from  this  figure.  This  indicates  that  the  material  is  not 
strain-rate  dependent,  although  it  exhibited  a  significant  amount 
of  hardening.  The  flow  stress  (5.5  kbars)  at  15  percent  strain  was 
more  than  twice  the  initial  yield  stress  (around  2  kbars).  The 
high  temperature  SHB  test  results  are  shown  in  Figure  7.  The 
stress- stra 1 n  curves  at  various  temperatures  show  the  reduction  in 
flow  stress  due  to  increase  in  temperatures. 
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TABLE  1 

SUMMARY  OF  SHB  TESTS 


i  — 

1 

TEST 

temperature 

1 

TYPE 

MATERIAL 

• 

STRAIN  RATE 

1  /  4  >  1 

wB-39o  1 

ROOM 

i 

TENSION 

1044  STEEL 

1 

1200 

HB-397  • 

ROOM 

\ 

TENSION 

1044  STEEL 

1 

1225 

H3-398 

ROOM 

1 

TENSION 

1044  STEEL 

1 

HB-399 

ROOM 

TENSION 

1044  STEEL 

700 

HB- 403 

ROOM 

1 

TENSION 

1044  STEEL 

1 

1  150 

HB-392 

ROOM 

1 

TENSION 

Cast  4140 

1 

HB-393 

ROOM 

1 

TENSION 

CAST  4140 

1 

H3-394 

ROOM 

1 

TENSION 

CAST  4140 

< 

H3-395 

ROOM 

1 

TENSION 

CAST  4140 

1 

\ 

HB-'OO 

ROOM 

« 

TENSION 

ARMCO  IRON 

l 

875 

1 

HB-401 

ROOM 

1 

TENSION 

ARMCO  IRON 

1 

250 

1 

HB-402 

ROOM 

t 

TENSION 

ARMCO  IRON 

1 

f— * 

O 

w 

1 

HB-404 

ROOM 

' 

TENSION 

ARMCO  IRON 

1 

1025 

1 

HB-405 

ROOM 

! 

TENSION 

ARMCO  IRON 

1 

800 

\ 

HB— 406 

ROOM 

1 

TENSION 

ARMCO  IRON 

1 

775 

1 

HB-407 

ROOM 

1 

TENSION 

ARMCO  IRON 

475 

1 

HB-406 

ROOM 

1 

TENSION 

ARMCO  IRON 

1 

4/5 

1 

HB-409 

1  ROOM 

1 

TENSION 

ARMCO  IRON 

1 

1300 

t 

HB-410 

!  ROOM 

t 

TENSION 

ARMCO  IRON 

\ 

1300 

1 

HB-41 1 

•  ROOM 

1 

TENSION 

ARMCO  IRON 

1 

1300 

1 

HB-413 

!  ROOM 

1 

COMP 

ARMCO  IRON 

1 

HB— 41 4 

!  ROOM 

1 

COMP 

ARMCO  IRON 

150 

1 

HB-41 5 

'  ROOM 

1 

COMP 

ARMCO  IRON 

1 

1600 

1 

HB-374 

!  Kuun 

■»  *■  •  t  nn 

ICN31UM 

w —  mo no 

rid  i  iu-v/u 

1 

540 

1 

1 

HB-375 

1  ROOM 

1 

TENSION 

Mar-M200 

1 

1020 
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TABLE  1  (continued) 


SUMMARY  OF  SHB  TESTS 

TEST 

TEMPERATURE  ' 

TYFE  1 

MATERIAL 

STRAIN  RATE  <  1  /  *  ) 

HB-376 

ROOM  ! 

TENSION  ' 

Mar-M200  * 

7S0 

HB-390 

ROOM  1 

TENSION  ' 

Mar-M300 

800 

HB-391 

ROOM  ' 

TENSION  1 

Mar-M300  ' 

700 

HB-296 

ROOM 

TENSION 

52125  STEEL' 

660 

HB-297 

ROOM  ' 

TENSION 

52125  STEEL' 

680 

HB -302 

700  C 

TENSION 

52125  STEEL1 

840 

HB-303 

620  C 

TENSION 

52125  STEEL' 

800 

HB-314 

ROOM 

TENSION 

52125  STEEL' 

525 

HB-315 

ROOM  ' 

TENSION 

52125  STEEL' 

57  5 

HB-316 

ROOM  ' 

TENSION 

52125  STEEL' 

1150 

■  •  ^ 

nD"o a  f 

ROOM  ! 

TENSION 

52125  STEEL ' 

1050 

HB-31B 

ROOM 

COMP 

52125  STEEL' 

200 

HB-319 

ROOM  ' 

COMP 

52125  STEEL' 

300 

HB-320 

ROOM  ! 

COMP 

52125  STEEL ! 

300 

HB -32 1 

ROOM  ! 

COMP. 

52125  STEEL1 

200 

HB-326 

ROOM  ! 

COMP 

52125  STEEL' 

400 

HB-327 

ROOM 

COMP 

52125  STEEL' 

400 

H8-334 

227  C 

TENSION 

52125  STEEL' 

725 

HQ-343 

445  C 

TENSION 

52125  STEEL1 

800 

HB-346 

700-750  C  ! 

TENSION 

52125  STEEL 

775 

HB-350 

HIGH  UNKNOWN  ! 

TENSION 

52125  STEEL' 

550 

HB-362 

ROOM  NECKING  ' 

TENSION 

52125  STEEL1 

HQ-363 

445  C 

TENSION 

52125  STEEL 

725 

HB-298 

ROOM 

TENSION 

NICKEL 

840 

HB-299 

ROOM  ' 

TENSION 

NICKEL 

840 
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TABLE  1  (continued) 
SUMMARY  OF  SHB  TESTS 


» - - “ - -  - - - - - - ------ - , - . - - - - —  I 


'  TEST 

TEMPERATURE 

1 

TYPE 

material 

1 

STRAIN  RATE  <  1  /s  )  1 

'  HB-300 

455  C 

1 

TENSION 

NICKEL 

l 

1120  ' 

•  HB-301 

960  C 

t 

TENSION 

NICKEL 

1 

1150  ' 

'  HB-322 

ROOM 

1 

COMP 

NICKEL 

1 

GOO  • 

'  HB-323 

ROOM 

\ 

COMP 

NICKEL 

1 

800 

'  HB-324 

ROOM 

1 

COMP 

NICKEL 

1 

1000 

'  HB-325 

ROOM 

t 

COMP 

NICKEL 

1 

1000 

!  HB-328 

ROOM 

1 

COMP 

NICKEL 

1 

BOO  1 

'  HB-329 

ROOM 

1 

COMP 

NICKEL 

t 

400 

!  HB-333 

!  210  C 

1 

TENSION 

NICKEL 

t 

775 

!  HB-342 

'  328  C 

* 

TENSION 

NICKEL 

1 

575 

•  HB-345 

1  713  C 

1 

TENSION 

NICKEL 

1 

850  ! 

!  HB-365 

■  252  C 

1 

TENSION 

NICKEL 

1 

525 

f  HB-366 

1  505  C 

1 

TENSION 

NICKEL 

1 

525  • 

•  HB-367 

!  675  C 

1 

TENSION 

NICKEL 

1 

575  • 

1  HB-36B 

1  244  C 

1 

TENSION 

NICKEL 

1 

1625  ! 

'  HB-369 

!  502  C 

1 

TENSION 

NICKEL 

1 

1625  ' 

!  H8-370 

1  700  C 

« 

TENSION 

NICKEL 

i 

1 600  ! 

'  HB -37 1 

•  ROOM 

1 

TENSION 

NICKEL 

1 

•  HB-372 

1  ROOM 

1 

TENSION 

NICKEL 

i 

!  HB-3BO 

'  ROOM 

1 

COMP 

VANADIUM 

■ 

100  ' 

!  HB-381 

'  ROOM 

1 

COMP. 

VANADIUM 

1 

120  1 

'  HB-382 

'  ROOM 

1 

COMP. 

VANADIUM 

1400 

'  HB-383 

'  ROOM 

1 

COMP 

VANADIUM 

1 

1300  1 

'  HB-366 

1  nuun 

1 

CQMP 

VANADIUM 

1 

1200 

!  HB— 389 

1 

'  ROOM 

i 

COMP 

VANADIUM 

1 

700  ! 
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TABLE  1  (concluded) 


SUMMARY  OF  SHB  TESTS 


TEST  - 

TEMPERATURE 

1 

TYPE 

1  MATERIAL 

•  STRAIN  RATE  < 1 /%  ) 

HB-335  ' 

ROOM 

1 

TENSION 

!  TUNGSTEN 

'  925 

HB-336  ' 

ROOM 

1 

TENSION 

'  tungsten 

!  675 

HB-337 

550 

C 

1 

TENSION 

!  TUNGSTEN 

1  450 

MB -338 

650 

C 

1 

TENSION 

'  TUNGSTEN 

'  500 

HB"339 

335 

C 

( 

TENSION 

'  TUNGSTEN 

*  400 

HB-3A7 

320 

c 

1 

TENSION 

!  TUNGSTEN 

1 

HB— 348 

195 

c 

1 

TENSION 

1  TUNGSTEN 

'  625 

HB-349 

95  C 

1 

TENSION 

'  TUNGSTEN 

!  500 

HB— 282 

120 

c 

» 

TENSION 

1  6061-T6 

•  850 

HB-283 

120 

c 

1 

TENSION 

•  6061-T6 

1 

HB— 284 

150 

c 

1 

TENSION 

!  6061-T6 

1  920 

HB-2BQ 

210 

c 

« 

TENSION 

'  6061-T6 

!  775 

HB-289 

El 

D 

TENSION 

!  6061 -T6 

'  750 

HB— 290 

345 

c 

I 

TENSION 

!  6061-T6 

!  800 

HB-291 

347 

c 

( 

TENSION 

!  6061 -T6 

!  800 

HB— 304 

190 

c 

I 

TENSION 

'  6061-T6 

•  880 

HB-305 

ROOM 

i 

COMP 

!  6061-T6 

'  800 

HB-353 

458 

c 

i 

TENSION 

!  606 1 -T6 

!  975 

HB-353 

458 

c 

i 

TENSION 

!  ALUMINUM 

!  975 

E  STRESS  Ckbo.0 


ARMCO  IRON 


STRAIN  RATE  (1/sec) 


Figure  4.  True  Stress  vs.  Log  r_  Plot  for  Armco  Iron 
at  Different  Strain  Levels. 
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TRUE  STRESS  (kb or) 


TRUE  STRESS  CKBAR) 


TRUE  STRAIN 


Figure  7.  True  Stress  vs.  True  Strain  Plots  at 
Elevated  Temperatures  lor  Nickel  200. 
(Results  are  for  Tension  SHB  Tests  at 
Average  Strain  Rates  of  600  to  1200  sec-^. ) 
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c .  W-2 _ Tungsten 

Both  room  and  high  temperature  SHB  tests  were 

conducted.  The  test  temperatures  ranged  from  room  to  650°C.  The 

elevated  temperature  results  were  obtained  at  strain  rates  of 
_  1  _  1 

400s  to  600s  .  The  tension  tests  at  room  temperature  were 

performed  at  strain  rates  which  ranged  between  250s  1  to  1250s  \ 

In  addition  to  these  SHB  tests,  two  quasi-static  tests  at  strain 
rates  of  0.001s  ^  and  1.0s  ^  were  conducted.  The  results,  shown  in 
Figure  8a,  indicate  that  the  flow  stresses  at  1.0s  are  higher 
than  at  0.001s  .  Thus,  the  flow  stresses  at  SHB  strain  rates  are 
expected  to  be  much  higher  than  under  quasi-static  tests  because  of 
the  rate  sensitivity  of  the  material.  However,  because  this 
material  failed  at  very  low  strains  (-0.04),  the  SHB  stress- strai n 
curves  could  not  be  resolved  due  to  the  mechanical  noises  as¬ 
sociated  with  the  SHB  test  setup  at  low  strain  values  as  can  be 
seen  in  Figure  8b.  Based  on  the  several  tests,  we  estimated  an 
approximate  stress  value  of  12-13  kbar  at  strain  rates  between 
500-1  000s'1  . 

The  high  temperature  test  results  are  shown  in 
Figure  9c.  The  flow  stress  dropped  significantly  at  elevated 
temperatures  while  the  strain  at  failure  increased  significantly. 

d .  Cl  008  Steel 

Wc  reported  the  quasi-static  and  SHB  room  tem¬ 
perature  results  earlier  in  Reference  1.  Several  high  temperature 
SHB  tensile  tests  were  conducted.  The  s tres3-strai n  curves  at  room 
temperature  and  at  670°C  are  shown  in  Figure  9d .  The  strain  rates 
in  these  tests  were  around  900s  1  . 

e .  60 6 1  -T6  Alumi num 

This  aluminum  is  a  strain  rate  independent 
material.  It  is  known  to  exhibit  slight  rate  dependency  above  a 
strain  rate  of  IQQQs  ^  .  The  SHB  and  quasi-static  test  results  for 
6061-T6  have  been  reported  earlier  in  Reference  10.  We  conducted 
SHB  tests  at  elevated  temperatures  using  the  induction  coil  heating 
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Figure  9.  True  Stress-True  Strain  Plots  at  Elevated  Temperatures 
for  (a)  6061-T6  Aluminum,  (b)  Nickel  200,  (c)~w-2 

Tungsten,  and  (d)  C1008  Steel. 
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technique.  The  results  are  shown  in  Figure  9a.  The  reduction  Jn 
strength  (flow  stress)  between  the  room  temperature  and  150°C  is 
not  significant.  Above  150°C,  the  flow  stress  reduced  rapidly  for 
increasing  temperatures.  This  feature  can  be  seen  from  the  earlier 
Figure  3. 

f .  Additional  Results 

We  conducted  necking  SHB  tests  on  1020  steel, 
606I-T6  aluminum,  HY100  steel,  C1008  steel,  W-2  tungsten,  and  52125 
steel.  In  Figure  10,  the  stress-strain  curves  beyond  the  onset  of 
necking  for  HY10C,  1020  steel,  and  606I-T6  aluminum  are  shown.  The 
data  for  other  materials  were  not  useful  due  to  preceding  failures 
of  the  specimen  before  any  significant  necking.  Detailed  descrip¬ 
tions  of  the  data  reduction  scheme  and  its  validity  can  be  found  in 
Reference  8.  Based  on  our  earlier  reported  results  for  these  three 
materials  [1],  HY100  steel  is  the  most  rate-dependent,  1020  steel 
is  less  rate  dependent,  and  6061-T6  aluminum  is  the  least  rate 
dependent  of  the  three.  Consideration  of  relative  rate  dependency 
aids  interpretation  of  the  stress-strain  curves.  In  Figure  10,  the 
3train  rate  is  indicated  at  various  points  along  the  curves.  The 
strain  rate  beyond  the  onset  of  necking  increased  continually  as 
the  neck- f ormat ion  progressed.  The  increase  in  flow  stress  beyond 
the  onset  of  necking  is  most  significant  in  the  most  strain-rate 
dependent  materials.  The  dotted  portions  of  the  curves  indicate 
the  estimated  strain- hardening  behavior  based  on  a  simple  power-law 
model  (0  -  Ken).  The  additional  stress  above  the  dotted  line 
indicates  the  level  of  strain- rate  sensitivity  in  the  metals. 

Since  606I-T6  aluminum  was  the  least  rate-dependent  material,  there 
i 3  relatively  little  increase  in  the  stress  due  to  strain  rate. 

The  additional  stress  in  1020  steel  was  higher  than  in  aluminum. 
HY100  was  the  most  rate  dependent,  and  the  increase  in  flow  stress 
is  the  most  significant.  When  the  strain  rate  increased  most 
rapidly,  the  flow  stress  also  showed  the  steepest  increase. 

We  also  acquired  SHB  tensile  and  compression  data 
for  several  other  metals.  The  stress-strain  curves  for  vanadium 
(compression  only),  52125  steel,  Mar-M200  steel,  Mar-M300  steel, 
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effective  stress  (kbar) 
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Figure  10. 


True  Stress-True  Strain  Plots  Obtained  from 
High  Speed  Photographs  of  Necking  Specimens. 


cast  4140  steel,  tantalum,  and  1044  steel  are  given  in  Figures  11- 
16.  No  reference  quasi-static  data  were  obtained  for  these 
materials . 

2.2  LONG  ROD  IMPACT  EXPERIMENTS 

The  need  for  developing  the  long  rod  experimental  technique 
in  which  a  flyer  plate  impacts  a  long  rod,  was  due  to  the  inade¬ 
quacy  of  SHB  test  technique  in  determining  the  initial  yield 
stress.  The  SHB  technique  is  very  useful  for  determining  flow 
stresses  at  plastic  strains  beyond  three  to  five  percent.  The 
plate-impact  technique,  on  the  other  hand,  is  very  useful  for  the 
determination  of  dynamic  yield  strength  (at  low  strains),  at  very 
high  strain  rate.  Thus,  it  is  very  difficult,  if  not  impossible, 
to  directly  compare  SHB  and  plate  impact  results,  especially  for 
rate-dependent  materials. 

2.2.1  Techni que 

A  new  technique  has  been  used  at  the  University  of 

Dayton  Research  Institute  which  is  based  on  the  long-rod  impact 

experiments  described  in  Kef.  11.  These  experiments  give  accurate 

values  for  the  yield  strength  of  the  specimen  at  strain  rates  of 
2-1 

above  10  s  ,  thus  bridging  the  gap  between  static  experiments  and 
the  plate-impact  experiments. 

A  schematic  representation  of  a  long-rod  experiment 
is  shown  in  Figure  17.  A  50  mm  diameter  steel  disk  is  used  for  a 
heavy  flyer  plate.  The  rod  is  12.7  mm  in  diameter.  The  gage  is 
inserted  at  a  distance  of  ten  diameters  away  from  the  impacted  end. 
Upon  impact,  both  elastic  and  plastic  waves  are  generated  in  the 
rod  and  propagate  towards  the  gage  station.  Various  calculations 
show  that  far  enough  from  the  impact  face  (some  ten  diameters  or 
so)  a  state  of  uniaxial-stress  exists  in  the  rod.  Thus,  at  the 
gage  location  we  expect  an  elastic  uniaxial-stress  wave  to  appear 
some  time  after  impact.  Thi3  time  is  determined  by  the  lon¬ 
gitudinal  sound  3peed,  Go  -  ✓  E/  p  ,  where  E  is  Young’s  modulus  and  p 
is  rod  density.  The  elastic  wave  amplitude  should  be  equal  to  the 
dynamic  yield  strength  of  the  rod  material.  The  rise  time,  which 
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Figure  12.  True  Stress-True  Strain  Plots  for  52125  Steel  under  Tension. 
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Figure  12.  True  Stress-True  Strain  Plots  for  52125  Steel  Under  Tension (concluded) . 


Stress-True  Strain  Plots  for  Marag: 


TRUE  STRESS 


1044  STEEL  W/ROLLING  ROOM  TEMP3 


Figure  14.  True  Stress-True  Strain  Plots  for  1044  Steel  Under 
Tension  Rolling  Direction. 
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Figure  15 


True  Stress-True  Strain  Plots  for  Pure  Tantalum 
Under  Tension. 
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True  Stress-True  Strain  Plots  for  Cast  4140  Steel 
at  k  =  1200s“l. 
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17.  A  Schematic  Diagram  of  Long  Rod  Experiment. 
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is  relatively  large  because  of  geometrical  dispersion,  can  be  used 
to  derive  the  (elastic)  strain  rate. 

The  fact  that  the  gages,  which  were  calibrated  under 
dynamic  uniaxial  strain,  can  be  used  to  measure  waves  in  rod3, 
which  experience  uniaxial  stress,  has  been  explained  in  Ref.  11. 

It  is  assumed  that  its  piezoresistance  behavior  is  the  same  as  in  a 
planar  shook  experiment.  The  two  parts  of  the  rod,  together  with 
the  gage  and  25  pm  thick  mylar  sheets  on  each  of  its  sides,  are 
glued  together  with  a  Hysol  0151  clear  epoxy,  in  a  special  fixture 
which  insures  good  axiality  of  the  parts  while  the  epoxy  is 
curing . 

For  a  rate- insens i ti ve  material,  such  as  2024-T3 
aluminum,  the  measured  elastic-wave  amplitude  is  exactly  the  same 

as  the  static  yield  strength  of  this  material  (see  Reference  11). 

-3  -  1 

Thus  the  yield  strength  is  the  same  at  strain  rates  of  10  s 
Moreover,  from  HEL  measurements  on  this  alloy,  the  dynamic  yield 

fi  -  1 

strength  at  10  s  is  al30  very  close  to  the  static  value,  confirm¬ 
ing  the  rate  insensitivity  of  this  material. 

A  typical  gage  record  from  a  long-rod  experiment  is 
shown  in  Figure  18.  We  first  note  the  arrival  time  of  the  stress 
wave  (some  20  ps)  which  is  due  to  geometric  dispersion.  A  small 
amplitude  oscillation  is  superimposed  on  the  signal  in  a  similar 
way  to  that  observed  in  the  earlier  work  on  aluminum  rods  [11]. 

Experiments  were  conducted  on  three  steels:  1020, 
C1008  and  HY1Q0.  The  amplitude  for  the  elastic  waves  for  the  three 
steel  specimens  is  given  in  Table  2  together  with  the  relevant  data 
concerning  these  shots..  In  order  to  compare  these  results  with  the 
static  and  plate  impact  experiments,  we  listed  the  various  yield 
strengths  for  the  three  steels  in  Table  3- 
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Figure  18.  A  Typical  Stress  Gauge  Record  from  a 
Long  Rod  Experiment. 
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TABLE  2 

RESULTS  FROM  LONG-ROD  EXPERIMENTS 


Shot  No. 

Rod  Material 

Impact 

Velocity  (m/s) 

Y ^  (GPa  )* 

925 

Cl  008 

1  96 

3.7 

926 

HY1  00 

201 

8.7 

897 

1  020 

94 

5.5 

898 

1  020 

1  64 

5.2 

_ 

Initial  yield  stress, 


TABLE  3 

THE  YIELD  STRENGTHS  (IN  GP  a)  FROM  DIFFERENT  TESTS 


Material 

Plate  Impact 
(106s-1 ) 

Long  Rod 
(  102-103s_1  ) 

Static  Tests 

Oo'V1  ) 

C1008 

0.65 

0.47 

0.25 

HY100 

1  .03 

0.87 

0.73 

1  020 

0.65 

0.53 

0.25 

The  1020  and  C1008  are  relatively  soft  steels  com¬ 
pared  to  the  HY100  steel  (static  yield  strength  of  7.3  kbar  in 
HY100  as  compared  to  2.5  kbar  for  the  other  two  steels)  [1]. 

2.2.3  Gage-Signal  Analysis 

Analysis  of  the  gage  record  shown  in  Figure  18  is 
presented  in  this  section.  A  small  precursor  arrives  at  the  gage 
location  at  a  time  which  corresponds  to  a  longitudinal  wave 
velocity  «  [(K  +  4/3G)/p]  (where  K  and  G  are  bulk  and  shear 
moduli,  respectively).  The  value  of  is  6  mm/ps  and  is  higher 
than  the  sound  3peed  Cn  =  4.6  mm/ps  (CQ  =-  /E/p ,  where  E  is  Young's 
modulus).  Their  existence  shows  that  the  longitudinal  elastic 
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waves  do  not  damp  completely  even  at  these  large  distances  from  the 
impact  end. 

The  wave  reflections  in  the  long  rod  can  be  explained 
with  the  aid  of  a  Lagrangian  x-t  diagram  as  shown  in  Figure  19. 

Upon  impact,  an  elastic  precursor  (compression  wave)  travels 
through  the  rod.  It  arrives  at  the  gage  location  at  point  'A*  and 
reflects  from  the  free  end  at  'B'  as  a  release  wave  towards  the 
gage  section  at  point  '  C*.  Further  towards  the  impact  end,  this 
release  wave  encounters  the  slow-moving  plastic  boundary  at  point 
'  D'.  It  reflects  back  at  this  boundary  towards  the  gage  section  as 
a  second  elastic  compression  wave.  By  observing  the  time  of  ar¬ 
rival  of  this  reflected  wave,  we  can  determine  the  distance  which 
the  plastic  front  travelled  from  the  impact  point  until  it  met  the 
elastic  release  wave.  Thus,  for  C1008  steel,  this  interaction 
occurred  33  mm  away  from  impact;  while  for  the  HY100  it  was  at  44 
mm.  This  difference  is  consistent  with  the  fact  that  in  the  SHB 
experiments,  HY100  experiences  more  work  hardening  than  C1008  (see 
Ref.  1).  The  results  for  HY100  steel  obtained  from  the  manganin 
stress  gages  compared  extremely  well  with  the  results  on  a  similar 
steel  reported  by  Kuscher  et  al .  in  Reference  12. 

As  mentioned  earlier,  the  SHB  is  not  accurate  enough 
in  the  low  strain  region  to  provide  a  reliable  value  for  the 
initial  yield  stress.  Thus,  the  long-rod  experiment  fills  the  gap 
between  the  static  and  shock  regimes  and  complements  the  data  for 
SHB  measurements  which  give  s t ress- strai n  curves  at  plastic  strains 
(beyond  3-5  percent).  To  accurately  determine  the  strain  rate  for 
the  long-rod  experiment,  we  used  the  measured  wave  amplitude  and 
the  Young's  modulus  of  steel  to  get  the  strain  carried  by  the  wave 
(in  the  range  of  0.25-0.4  percent)  and  th<-  rise  time  of  the  wave 
(4-5  microseconds).  These  result  in  strain  rates  in  the  500~800s  1 
region.  The  main  contribution  of  the  long-rod  experiment  is  that 
with  a  relatively  simple  and  straightforward  technique  we  can 

2 

measure  the  yield  strength  of  the  specimen  at  strain  rates  of  10 

3  - 1 

to  10  s  .  These  values  can  then  be  used  to  complement  SHB  data 
which  lack  accuracy  in  the  low  strain  region. 
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2.3  EXPANDING  RING  EXPERIMENTS 


The  expanding  ring  technique  presents  an  attractive  supple¬ 
ment  to  conventional  tensile  test  techniques.  Relatively  high 

14  —  -j 

rates  can  be  achieved,  >10  s  .  However,  the  most  attractive 
feature  of  ring  tests  is  that  the  data  analysis  is  free  from  as¬ 
sumptions  involving  wave  propagation. 

In  the  expanding  ring  test,  an  axisymmetric  thin  short  ring 
is  driven  racially  outward.  The  ring  may  be  launched  explosively, 
with  a  magnetic  field  or  with  an  exploding  wire  via  a  driver 
cylinder.  Most  past  work  with  expanding  rings  has  been  designed  to 
measure  ultimate  strain  at  fracture.  The  ultimate  strain  is  given 
by: 


ef  =  in  rf/rQ  (1) 

Here  rQ  and  are  the  initial  radius  of  the  cylinder  and  the 
radius  of  failure,  respectively.  When  the  average  velocity,  r,  is 
monitored,  strain  rate  may  be  calculated  from; 

e  =  r/ r  ( 2 ) 

The  strain  in  the  ring  can  be  determined  from  the  instan¬ 
taneous  diameter.  The  condition  that  the  ring  is  in  free  flight 
permits  the  stress  in  the  ring,  ,  to  be  determined  by  measuring 
the  radial  deceleration: 


oe  =  -prr  (3)  ' 

The  experimental  setup  to  expand  the  ring  dynamically  con¬ 
sists  of  a  capacitor  bank  with  detonating  wire  assembly,  basically 
as  described  by  Rajendran  and  Fyfe  [13].  An  average  deceleration 
can  be  measured  from  the  record  using  a  streak  camera  -  flaahlamp 
optical  setup.  We  introduced  several  modifications  to  both  the 
mechanical  operation  and  the  electronic  circuitry.  This  includes: 
(a)  modification  of  the  rail  spark  gap  reducing  its  separation  to 
1/4"  and  hookup  of  its  trigger  electrode  by  a  voltage  divider,  (b) 


replacement  of  the  kv  trigger  module,  (c)  design  and  installation 
of  a  safety  switch,  which  allows  testing  of  pref ire/no  fire  without 
discharging  the  hank,  and  (d)  shadow  graph  recording  system,  using 
a  Xenon  flashlight  and  a  streak  camera.  Tests  conducted  showed  a 
complete  detonation  of  a  copper  wire  0.85  mm  thick,  discharging  19 
kv  from  the  87  pf  capacitor  bank  at  an  energy  equivalent  of  about 
1  6  KJ . 

Currently,  the  expanding  ring  experimental  setup  is  being 
routinely  used.  We  conducted  several  tests  on  6061 -T6  aluminum  and 
1100-0  aluminum  to  demonstrate  the  capabilities  of  the  system.  In 
Figure  20,  a  sample  streak  record  is  shown.  On  this  Air  Force 
contract,  we  updated  our  expanding  ring  facility  to  acquire  high 
strain  rate  properties  routinely.  Future  efforts  will  include 
acquiring  data  for  several  metals. 

2.4  PLATE  IMPACT  EXPERIMENTS 

In  the  plate  impact  test,  a  flat  flyer  plate  is  impacted 
against  a  target  plate  at  a  high  velocity.  The  flyer  and  target 
may  be  of  the  same  or  different  material.  Compressive  stresses  are 
produced  and  transmitted  rapidly  from  the  plane  of  impact  to  the 
adjacent  stress  free  areas  of  the  material  in  the  form  of  a  stress 
pulse.  Many  discussions  of  planar  impact  loading  are  available 
[14,15]. 

Plate  impact  test  provide  a  loading  path  that  is  very  dif¬ 
ferent  from  conventional  SHB  compression  or  tension  tests.  The 
deformation  is  that  of  one-dimensional  strain,  and  the  mean  stress 

is  generally  very  high  compared  to  that  of  SHB  tests.  Strain  rates 

5  - 1 

are  10  a  or  higner.  The  material  undergoes  compression  im¬ 
mediately  followed  by  tension  due  to  a  reflected  wave  from  the  rear 
surface  of  the  target.  Plate  impact  experiments  are  essential  for 
calibrating  and  validating  high  strain  rate  material  models  that 
aspire  to  general  applicability.  Specifically,  plate  impact  data 
may  be  interpreted  to  infer  compression  and  tensile  yield  strengths 
and  failure  parameters. 
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The  plate  impact  tests  were  conducted  with  three  objectives: 
(1)  determination  of  Hugoniot  Elastic  Limit,  (2)  determination  of 
the  unloading  path  from  the  free  surface  velocity  history,  and  (3) 
determination  of  the  threshold  conditions  associated  with  onset  of 
spall  fracture. 


Impact  induces  an  elastic  shock  and  a  plastic  shock  in  the 

target.  The  amplitude  of  the  elastic  shock  is  oHEL.  The  Hugoniot 

Elastic  Limit,  oHEE,  is  the  maximum  stress  in  the  elastic  regime 

for  one-dimensional  wave  propagation.  This  stress  level  is 

generally  treated  as  a  material  property,  and  above  this  level  the 

material  flows  plastically.  The  stress  o.,,-,,  can  be  determined  from 

HRL 

the  experimentally  obtained  free  surface  velocity  of  the  target 
that  corresponds  to  the  elastic  shock  velocity,  uHEE 


°HEL 


1/2  PCL 


UHEL 


(4) 


where  is  the  elastic  sound  speed  in  one-dimensional  strain.  The 
high  strain  rate  yield  strength  XQ  can  then  be  calculated  from  the 
relationship: 


HEL 

0  =  (Ji  ♦  1) 

^  2G  3  1 


(5) 


where  G  is  the  3hear  modulus  and  K  is  the  bulk  modulus. 

2.4.1  Results 

Plate  impact  tests  were  conducted  on  W-2  tungsten, 
1215  steel,  5083-0  aluminum,  and  Armco  Iron.  In  addition,  for 
precursor  decay  analysis,  several  additional  plate  impact  tests  on 


C1008  steel  were  conducted.  Results  in  terms  of  o 
stress  for  these  metals  are  given  in  Table  4. 
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A  number  of  new  or  improved  techniques  for  high  strain  rate 
measurements  were  developed.  These  are  described  in  detail  in  this 
section.  Various  applications  are  also  discussed. 
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TABLE  4 


OBSERVED  HEL  AND  SPALL  THRESHOLDS  FOR  CRACK  FORMATION 


Material 

HEL 

(Kbar) * 

Spall  Threshold 
Impact  Velocity 
(M/S) 

Spall 

Stress,  o 
(kbar)  s 

1215 

Steel 

14.4+02 

230 

8.0  ^ 

W2 

Tungsten 

24.0 

210 

7.4 

5083-0 

Aluminum 

2.9 

426 

16.0 

1  Kbar  =  14504  psi 

=  100  MPA 

9  2 

=  10  dynes/cm 


Incipient  spall;  low  spall  strength  because  spall  signal 
was  due  to  precursor  voids  in  the  target. 


2.5.1 


Manganin  G auge  Technique 


P iezoresistive  gauges  have  been  used  extensively 
during  the  past  20  years  as  in-material  stress  transducers  in 
planar  shock  wave  studies.  Manganin  and  Ytterbium  are  the  most 
popular  among  these  transducers.  Many  workers  determined  calibra¬ 
tion  curves  for  manganin,  in  the  0-200  kb  range,  in  terms  of 
relative  resistance  change  versus  shock  stress.  These  are  reviewed 
in  References  2,  16,  ana  17-  There  are  two  important  difficulties 
with  the  use  of  piezores.stance  stress  gauges  in  dynamic 
experiments.  First,  there  is  the  interpretation  of  data  for  load¬ 
ing  conditions  whi.  ;h  differ  from  the  standard  calibration 
experiments.  The  most  prominent  example  is  the  use  of  piezoresis- 
tive  gauges  to  measure  lateral  stress  in  i-D  experiments.  A  second 
difficulty  arises  from  the  appearance  of  a  substantial  resistive 
hysteresis  in  the  gauge  response  upon  unloading.  For  example,  it 
can  be  shown  that  wire  gauges  should  have  a  different  response  than 
foil  gauges. 

We  used  and  refined  a  rather  simple  analytic  elasto- 
plastic  model  for  the  gauge  response  which  is  a  direct  continuation 
of  the  work  [2].  The  present  approach  is  general  and  can  be  used 
for  any  piezoresistance  gauge.  However,  manganin  gauges  are  the 
only  gauges  for  which  we  have  a  large  amount  of  data  so  that  most 
of  this  section  deals  with  these  gauges.  Specifically  we  deal  with 
three  issues:  (1)  accounting  for  the  differences  between  the 
various  calibration  curves  for  manganin  gauges,  (2)  demonstrating 
the  importance  of  the  yield  strength  of  the  gauge  material  In 
determining  its  resistive  hysteresis  upon  unloading  and  accounting 
for  the  higher  response  of  annealed  gauges,  and  (3)  the  response  of 
transverse  foil  gauges  which  measure  the  lateral  stresses  in  plane 
shock  wave  experiments. 

The  well-known  logarithmic  derivative  of  the  electri¬ 
cal  resistance  of  a  conducting  element  having  length  i,  cross 
section  A,  and  resistivity  p  is  given  by: 

R  =  p|  (6) 
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where  e^,  are  longitudinal  (in  the  current  direction)  and  cross- 
sectional  strains,  defined  positive  in  compression.  For  the 
various  configurations  which  one  encounters  in  1-D  experiments,  the 
gauge  is  emplaced  so  that  its  length,  width,  and  thickness  coincide 
with  the  three  principal  directions  of  the  stress  tensor  which 
defines  the  shocked  state.  This  point  is  very  important  because  it 
means  that  the  following  analysis  is  not  general  enough  to  include 
cases  where  the  gauge  i3  emplaced  in  an  arbitrary  orientation 
relative  to  the  shock  front.  Thus,  in  conventional  plate  impact 
teat,  shear  stresses  do  not  act  on  the  gauge  element. 

The  relative  resistivity  change  can,  in  principle, 
depending  on  the  plastic  work  done  on  the  gauge  through  the  lattice 
defects  which  are  generated  by  the  shock  wave.  These  extra  defects 
were  considered,  for  some  time,  to  be  the  main  reason  for  the 
resistive  hysteresis  in  manganin  gauges.  Recently,  Chen  et  al. 

[18]  found  that  for  manganin  these  defects  do  not  contribute  sig¬ 
nificantly  to  the  residual  resistivity.  This  confirms  the 
assumption  in  Reference  19  about  the  "mechanical"  nature  of  the 
hysteresis  in  manganin  rather  than  the  extra  defect  explanation. 

For  such  a  gauge  element,  which  is  aligned  along  the  z  direction, 
the  general  expression  (see  References  20  and  21)  for  the  resis¬ 
tivity  change,  is  given  by 
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where  ir  are  the  well-known  piezoresistance  coefficients  and  o ^ 
are  the  principal  stresses.  The  term  does  not  appear  in 

Equation  (7)  because  no  shear  stresses  act  on  the  gauge  element. 
Figure  21  shows  the  position  of  the  gauge  element  in  relation  to 
the  shock  front  and  the  principal  axes. 
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We  next  decompose  o_  to  its  hydrostatic  and 
deviatoric  terms  (o  =  p  +  o’)  and  the  following  relations  can  be 

Z  Zt 

obtained: 

=  (l711  "  1T12)  (P  +  °z)  +  3*12  P 

“  +  2  it  ^  2 )  P  +  ~  n12)  °2  (d) 

■  ”HS  •  p  *  11  shear  °i  ‘  l^'HS  *  se'f 

where  e^  is  the  elastic  deviatoric  strain  along  the  current  direc¬ 
tion  (z)  and  s  is  an  empirical  constant.  The  first  term  on  the 
right-hand  side  depends  only  on  the  volumetric  strain  of  the  gauge 
ey  and  can  be  determined  from  hydrostatic  compression.  The  second 
term  is  a  deviatoric  contribution  which  is  proportional  to  the 
deviatoric  part  of  the  elastic  strain  along  the  gauge  length  . 

We  now  make  the  usual  assumption  that  if  the  gauge 
element  is  much  longer  than  both  its  width  and  thickness,  then  the 
strain,  e^,  along  the  gauge  length  is  zero.  Since  cy  -  +  eA,  we 

get  for  the  relative  resistance  change  is: 
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This  equation  is  the  general  result  which  shows  that 
the  response  of  piezoresisti ve  gauges  is  mairly  determined  by  the 
hyrdodynamic  pressure  through  the  (Ap/pQ)HS  and  ev  terms>  The 
third  term  (s  •  e^)  is  a  deviatoric  one  and  its  magnitude  depends 
on  the  yield  strength  of  the  gauge. 

In  Reference  1,  the  method  to  obtain  s  for  manganin 
by  using  the  measured  relative  resistance  changes  under  both 
hydrostatic  and  uniaxial  strain  conditions  has  been  shown.  We  used 
the  compression  curve  for  manganin,  as  given  in  Reference  22: 

p  -  1 1 60 e v  +  41 2 0 e y ,  (10) 
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where  p  is  the  kbar.  The  most  accepted  values  for  the  piezoresis- 

-3 

tance  coefficient  under  hydrostatic  conditions  is  2.4  x  10 
(Q/£l)/kb.  The  1-D  coefficient  found  in  Reference  23  for  the  elas- 

_3 

tic  range  of  gauge  response  is  1.95  x  10  (£l/ft)/kb.  These  values 

result  in  a  value  of  s  -  0.18  for  the  empirical  constant  in  Equa¬ 
tion  (9).  Using  this  value  together  with  the  compression  curve, 
(Equation  (10)),  we  finally  get  for  manganin: 

(jp)  -  3.45  (1  ♦  2.8?ey)  +  o.lSe^  (11) 

0  eL  *  0 

The  stress-time  history  can  be  obtained  using  the  above  calibration 
relationship . 

The  analysis  presented  above  demonstrates  the  impor¬ 
tance  of  the  yield  strength  of  the  gauge  material  in  determining 
its  resistance  change  for  given  stress  and  strain  conditions.  We 
can  clearly  see  how  a  high  value  of  ¥  and  especially  one  which 
increases  with  pressure,  can  affect  the  value  of  hydrodynamic 
pressure  (p)  in  the  gauge  element,  which  is  the  main  factor  in  the 
gauge  response.  The  importance  of  the  yield  strength  can  be  fur¬ 
ther  demonstrated  by  considering  its  influence  on:  (1)  the 
resistive  hysteresis  of  the  gauge  and  (2)  the  response  of  annealed 
gauges.  In  the  next  subsections  we  elaborate  on  these  issues.  By 
observing  these  issues,  in  light  of  the  elasto-plastlc  response  of 
the  gauge  material,  we  get  a  better  understanding  of  the  dynamic 
loading  process.  These  points  are  also  discussed  in  Reference  2. 

2.5.2  Transverse  Gauge  Technique 

One  of  the  most  important  issues  in  the  study  of 
dynamic  behavior  of  solids  under  shock  loading  is  the  determination 
of  the  lateral  stress  in  the  shocked  specimen.  The  only  direct 
technique  for  such  a  measurement  uses  thin  piezore3istance  stress 
gauges  which  are  inserted  laterally  in  the  shocked  target  as  shown 
schematically  in  Figure  22. 

In  Reference  2  oast  work  on  direct  measurement  of 
lateral  stress  with  piezoresistive  gauges  has  been  reviewed.  The 


47 


SHOCK  DIRECTION-X 


Figure  22  -  The  Principal  Stresses  for  a  Transverse 
Stress  Gauge. 
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approach  used  in  this  program  is  baaed  on  that  work.  Our  technique 
is  also  described  in  Reference  24. 


In  order  to  develop  an  analytical  calibration  for  the 
transverse  gauge,  the  stress  across  the  thickness  of  the  gauge  is 
assumed  to  be  equal  to  the  lateral  stress  in  the  target  -  oy.  This 
assumption  is  actually  the  basic  one  behind  the  lateral  stress 
measurement  technique  with  piezoresisti ve  gauges,  and  it  is  jus¬ 
tified  for  very  thin  (foil)  gauges.  It  is  also  assumed  that  the 
strain  in  the  gauge  along  the  shock  direction  is  the  same  as  the 
uniaxial  strain  in  the  surrounding  specimen  -  ex>  This  means  that 
the  width  of  the  gauge  is  decreased  by  ex  (the  gauge  is  bonded 
perfectly  to  the  specimen).  To  justify  this  perfect  bonding  as¬ 
sumption,  it  is  best  to  refer  to  a  related  work  in  which  we  have 
shown  that  a  thin  wire,  or  foil,  which  is  placed  longitudinally  in 
the  target  (along  3hock  propagation  direction)  undergoes  the  same 
3train  as  its  surroundings. 

The  next  step  is  to  realize  that  the  elastic  con¬ 
stants  of  the  gauge  and  matrix  can  be  very  different  (as  in  the 
case  of  manganin  in  PMMA) .  Thus  the  assumption  of  equal  ex  means 
that  the  stress  in  the  gauge  along  its  width  (ox)  may  be  very 
large.  This  can  easily  be  seen  if  we  write  the  elastic  expression 
for  ox: 


ox  -  (K  -  |g)  ey  +  2Gex 


(12) 


where  ey  is  the  volumetric  strain  of  the  gauge  and  K  and  G  are  its 
bulk  and  3hear  moduli,  respectively.  Thus,  the  relatively  large 
strains  in  shock-loaded  PMMA  (about  1  percent  per  kbar)  will  result 
in  very  large  stresses  in  the  manganin  gauge  because  of  the  2Gex 
term  (2G  -  870  kbar;  see  Reference  22).  This  means  that  the  elas¬ 
tic  behavior  of  the  transverse  gauge  can  be  confined  to  a  very 
small  stress  range.  For  manganin  in  PMMA  this  range  is  limited  to 
shock  stresses  of  leas  than  about  0.5  kbar.  The  values  that  oy  can 
take  are  then  determined  by  the  value  of  oy  and  the  yield  criterion 
of  the  gauge  material.  Assuming  a  Von-Misses  (or  Tresca)  yield 
criterion,  we  thu3  write  for  o  in  the  gauge: 
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where  Y  is  the  yield  strength  of  the  gauge  material. 

There  remains  the  question  of  the  magnitude  of  the 
third  stress  component  along  the  gauge  length  -  o„.  Since  this  is 
the  largest  dimension  of  the  gauge  we  may  assume,  as  do  others, 
that  the  strain  along  this  direction  is  zero.  As  for  the  stress 
along  the  z-axis  we  assume  that  it  is  equal  to  oy  in  the  gauge. 

This  assumption  can  have  some  justification  if  we  consider  the 
calculations  of  Gupta  [25]  which  show  that  the  two  lateral  stresses 
in  the  gauge  ( oy  and  oz)  are  nearly  equal. 

Thus  the  stress  and  strain  tensors  which  describe  the 
transverse  gauge  are: 


We  note  here  that  because  ex  can  be  much  larger  than  ey  (as  in  the 
case  of  manganin  in  PMMA)  the  value  of  the  3train  through  the 
thickness  of  the  gauge  (ey)  can  be  negative.  Thus  the  thickness  of 
the  transverse  gauge  can  increase  under  shock  compression.  This  is 
also  what  Gupta  [25]  found  in  his  calculations  for  Ytterbium  in 
PMMA . 

After  having  the  principal  stress  in  the  gauge,  we 
can  write  an  expression  for  the  hydrostatic  pressure  which  prevails 
in  it: 


( 1  5) 
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The  next  step  is  to  consider  the  calibration  curve  of 
the  piezoresisti ve  gauge  in  terms  of  its  relative  resistance  change 
versus  its  volumetric  strain  (ev)  as  in  Equation  (11)  for  manganin. 
This  curve  can  also  be  given  by  the  hydrostatic  pressure  in  the 
gauge  with  the  aid  of  the  dynamic  compression  curve  of  the  gauge 
P  -  f ( ev )  as  in  Equation  (10).  Thus: 

=  f(ey)  «  g(P),  (16) 

where  AR/Rq  is  the  relative  resistance  change  of  the  gauge. 

For  manganin  it  was  found  (see  Reference  21)  that 
Equation  (16)  takes  the  form  given  by  Equation  (11).  Thu3  ev  is 
given  by  solving  this  quadratic  relation: 


ey  =  -0.1774  +  0 . 03042  *  j-gj  (17) 

where  A  AR/Rq. 

The  data  reduction  process  involves  only  two  steps. 
First,  use  the  measured  A  value  in  Equation  (17)  to  obtain  the 
value  of  ey.  Then  use  ey  in  Equation  (10)  to  obtain  the  value  of 
the  pressure  in  the  gauge  (P).  The  second  step  is  to  use  this 
value  of  P  in  Equation  (15)  to  get  the  value  of  Oy.  In  order  to 
allow  for  the  changes  of  Y  with  peak  pressure  we  use  the  Y  -  Y(P) 
curve  in  Reference  26. 

2.5.3  High  Pressure  Measurements 

The  direct  determination  of  rarefaction  wave 
velocities  in  shock-loaded  solids  is  important  for  high-pressure 
equation  of  state  data.  Various  experimental  techniques  have  been 
attempted  in  the  past  20  years,  and  these  are  summarized  in 
Reference  27.  Currently,  the  most  widely  used  technique  employs 
optical  flash  gaps.  This  technique  has  since  been  used  to  measure 
longitudinal  wave  velocities  for  a  variety  of  shocked  materials. 

In  this  technique,  a  relatively  thin  flyer  plate  strikes  a  stepped 
target  plate.  The  time  duration  (At)  of  the  shock  level  at  each 
step  is  recorded.  This  is  done  by  backing  the  stepped  target  with 
transparent  materials  which  emit  large  amounts  of  thermal  radiation 


(when  shocked)  and  using  optical  analyzers  to  record  these  radia¬ 
tion  versus  time  histories.  From  the  different  time  durations 
versus  the  position  of  the  optical  analyzer  (step  thickness)  one 
then  determines  the  catch-up  ratio  (R).  This  is  the  ratio  of 
target  to  flyer  thickness  where  the  release  wave  from  the  flyer’s 
back  face  overtakes  the  shock  wave  in  the  target.  The  longitudinal 
sound  velocity  (the  leading  edge  of  the  rarefaction  fan)  is  given 
by 
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(18) 


where  pQ  and  p  are  initial  and  shock-compressed  densities  of  the 
specimen  and  U3  is  its  shock  velocity. 

The  optical  technique  has  several  disadvantages  for 
use  at  the  University.  It  requires  extremely  sensitive  or  fast 
optical  recording  apparatus,  intricate  target  fabrication,  and  use 
of  noxious  gasses.  Therefore,  a  new  technique  was  developed  using 
manganin  gauges  as  the  time-of-arri val  sensors  (Reference  28).  The 
experimental  configuration  is  shown  schematically  in  Figure  23.  A 
thin  flyer  plate  impacts  a  stepped  target  plate.  At  the  center  of 
each  step  is  placed  a  thin  (5  ym  thick)  manganin  stress  gauge 
backed  by  a  thick  epoxy  disk.  For  metallic  targets,  a  thin  Mylar 
sheet  (10-25  ym)  is  placed  between  the  gauge  and  the  target.  All 
these  elements  (Mylar  sheet,  gauge,  and  epoxy  disk)  are  then  bonded 
together  and  to  the  target  with  epoxy  under  a  press  in  order  to 
minimize  bond  thickness  and  air  bubbles.  Thus,  the  gauge  is  sur¬ 
rounded  by  a  uniform  material  (epoxy)  and  is  close  enough  to  the 
interface  to  measure  the  stres3-time  history  there.  By  embedding 
the  gauge  in  plastic,  we  greatly  improve  its  time  resolution,  which 
is  determined  by  the  gauge  thickness  and  is  of  the  order  of  a  f ew 
nanoseconds . 

There  is  another  advantage  to  the  backface  over  the 
in-material  configuration  which  is  very  important  for  the  present 
application.  The  embedded  gauge  is  limited  to  shock  stresses  of 
200  kb  or  so.  This  is  so  because  most  polymers  (epoxy  and  Mylar 


52 


Figure  23.  Schematic  Representation  of  the 
Impact  Configuration. 
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included)  become  electrically  conducting  at  these  pressures  and 
they  short  circuit  the  gauge  and  its  leads.  With  the  backface 
configuration  the  stresses  in  the  specimen  can  be  very  high  (of  the 
order  of  megabar)  before  the  stresses  in  the  epoxy  backing  reach 
200  kb.  This  is  especially  true  for  materials  with  high  acoustic 
impedances  (like  copper,  steels,  tungsten,  etc.). 


Thus,  the  backface  configuration  for  manganin  gauges 
is  most  appropriate  for  catch-up  ratio  determination  because  it 
overcomes  the  two  difficulties  of  the  in-material  gauge  configura¬ 
tion  (long  rise  times  and  limited  stress  range).  Manganin  gauges 
have  two  advantages  over  the  optical  analyzers  as  stress-time 
sensors.  First  of  all,  they  are  much  simpler  to  use  and  their 
experimental  assembly  involves  fewer  components  and  is  less  time 
consuming.  Second,  they  are  calibrated  for  the  entire  stress  range 
of  0-200  kb  both  under  loading  and  unloading  conditions.  This 
enables  us  to  identify  the  arrival  of  other  release  waves  with  much 
more  certainty  than  with  the  optical  analyzers.  Moreover,  from  the 
measured  stress-time  history  at  the  specimen-epoxy  interface  we  can 
determine  the  spall  strength  of  the  target. 


In  order  to  demonstrate  the  usefulness  of  manganin 
gauges  as  time-of-arrival  sensors  we  performed  a  few  symmetric 
plate  impact  experiments  with  0FHC  copper  flyers  and  targets. 

Gauges  were  used  as  described  in  Section  2.5.1.  Their  main  use  in 
the  present  study  is  as  time-of-arrival  gauges.  Details  of  these 
experiments  can  be  found  in  Reference  28. 

Figure  24  3hows  the  extrapolation  of  the  At  vs  gauge 

location  points  to  obtain  the  catch-up  distance  for  this  shot.  The 

impact  velocity  was  1042  m/s.  The  value  of  15.8  mm  corresponds  to 

a  catch-up  ratio  of  R  -  6.15  (flyer  thickness  was  2.57  mm).  This 

value  together  with  the  measured  shock  velocity  of  U  -  4.65  mm/ys 

s 

are  then  inserted  in  Equation  (18)  to  give  a  value  of  -  5.73 
mm/ys  for  the  longitudinal  sound  speed  of  copper  compressed  to  220 
kb.  This  value  is  in  excellent  agreement  with  the  value  of  5.7  mm/ys 
obtained  by  Morris  et  al .  [293.  The  details  of  the  tests  on  1215 
steel  are  summarized  in  Table  5. 
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Linear  Extrapolation  to  Find  the 
Catch-up  Distance. 


T ABLE  5 


SUMMARY  Ob’  THE  RESULTS  ON  1215  STEEL 
A.  IN  MATERIAL  GAUGE  CONFIGURATION;  IMPACTOR  -1215  STEEL 


Shot  No. 

Impact or 
Thickness 
(mm) 

Impactor 
Velocity 
( m/s ) 

Target 

PMMA 

(mm) 

Shock 
Stress 
(kbar ) 

7-1027 

10.2 

31  3 

6+6 

51  ±2 

7-1028 

10.2 

446 

6+6 

87  ±3 

7-1039 

9.  52 

71  9 

6.35+6.35 

1  36±3 

7-1040 

9.52 

870 

6.35+6 . 35 

1  64±3 

7-1041 

9.52 

1  21  6 

6.35+6.35 

222  ±5 

B. 

BACK  GAUGE  CONFIGURATION;  IMPACTOR  -  1215  STEEL 

7-1015 

2.54 

231 

4  +  8 

4  6±2 

7-1029 

2.54 

901 

4  +  8 

1  60±5 

7-1000 

2.54 

101  3 

4  +  8 

1  7  9±5 

7-1006 

2.54 

1  326 

4+8 

238*5 

7-1020 

4.0 

1  6  6 1 

4  +  8 

304±5 

7-1051 

4.0 

1000 

8.0 

1  7  7  ±  5 

7-1052 

4.0 

1  4  90 

8.0 
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2.5.4  Shear  Wave  Measurements 


Detection  of  shear  waves  is  one  of  the  more  difficult 
problems  in  shock  wave  diagnostics.  During  this  program,  we 
developed  a  shear  wave  generation  technique,  and  showed  that  man- 
ganin  gauges  could  be  used  to  detect  shear  wave  arrivals  [30]. 

A  piezoresistance  gauge  was  used.  It  was  calibrated 
under  both  pressure  and  shear  conditions,  in  order  to  have  a  quan¬ 
titative  measurement  of  both  the  compressional  and  shear  wave 
amplitudes.  In  addition,  the  time  of  arrival  of  the  shear  wave  at 
the  gauge  location  can  be  used  to  measure  the  pressure  dependence 
of  the  shear  wave  velocity,  since  this  wave  follows  the  faster 
compressional  wave.  One  should  emphasize  that  at  present,  the  use 
of  piezoresistance  gauges  in  shook  wave  studies  is  limited  to  the 
more  simple  one-dimensional  configuration.  In  this  configuration, 
there  are  no  shear  stresses  acting  on  the  gauge  and  the  only 
relevant  component  of  the  stress  tensor  is  that  along  the  shock 
propagation  direction.  This  stress  component-  in  the  gauge  element 
is  equal  to  that  in  the  target,  thus,  enabling  the  measurement  of 
shock  stresses.  In  order  to  use  a  piezoresist ive  element  of  the 
measurement  of  a  shear  stress  component,  one  has  to  use  a  gauge 
which  has  been  calibrated  under  pure  shear  conditions.  We 
demonstrated  the  feasibility  of  this  technique  by  an  oblique-impact 
experiment  with  a  rnanganin  gauge  embedded  in  an  alumina  target. 

The  experimental  assembly  is  shown  In  Figure  25. 

The  piezoresistance  element  should  ideally  have  a 
large  and  well-calibrated  shear  coefficient.  Folycrystall ine 
solids  have  a  single  direct  shear  piezoresi stant  coefficient  which 
is  term  The  gauge  should  also  have  a  well-calibrated  pressure 

coefficient  in  order  to  measure  the  compressional  wave  amplitude 
directly.  The  proper  gauge  should  be  embedded  between  two  specimen 
disks  in  the  same  way  as  regular  piezoresistance  gauges  are  in¬ 
stalled  in  plate-impact  experiments.  Gupta's  technique  [25]  to 
embed  the  gauge  in  very  small  grooves  in  the  target,  without  using 
epoxy  and  relying  only  on  the  friction  between  the  two  disks,  might 
prove  advantageous  here. 
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GAUGE  LEADS 


Figure  25.  A  Schematic  Representation  of  an  Oblique 
Impact  Experiment  with  a  Piezoresistance 
Gage . 
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When  the  oblique  impactor  hits  the  target,  it  gener¬ 
ates  two  waves  which  propagate  away  from  the  impact  face.  The 
first  is  the  usual  uniaxial-strain  compress ional  wave.  It  will  be 
recorded  by  the  gauge  as  described  in  Section  2.5.1.  The  slower 
shear  wave,  moving  in  an  already  compressed  material  will  cause  an 
additional  resistance  change  for  the  gauge  at  some  later  time 
(depending  on  the  shear  wave  velocity).  The  resistance  change  will 
either  increase  or  decrease  the  current  resistance  of  the  gauge, 
depending  on  the  sign  of  n^.  If  is  large  enough,  this  addi¬ 

tional  resistance  change  can  be  used  to  accurately  measure  the 
amplitude  of  the  shear  wave.  The  arrival  time  of  the  shear  wave 
can  be  used  to  measure  the  pressure  dependence  of  its  velocity. 

This  can  then  be  used  to  determine  the  pressure  dependence  of  the 
shear  modulus  of  the  target  material  which  i3  an  important 
parameter  for  various  applications. 

The  advantage  of  the  proposed  technique  is  that  with 
a  single  gauge,  both  the  compressional  and  shear  wave  profile  can 
be  directly  determined.  Thus,  the  effects  of  shear-related  dynamic 
processes,  like  phase  transformation  and  spall,  can  be  studied  more 
closely  by  varying  the  amount  of  shear  continuously.  A  drawback  to 
the  proposed  technique  is  that  at  the  present  time  (to  our 
knowledge)  no  piezoresistance  gauge  commonly  used  in  shock  experi¬ 
ments,  has  been  calibrated  for  shear  response.  Our  guess  is  that 
of  the  two  most  commonly  used  gauges,  manganin  and  ytterbium,  the 
latter  will  be  more  adequate  for  shear  measurements  because  it  has 
much  higher-pressure  coefficient  (ir^  and  tt12).  The  only  shear 
gauges  we  know  of  are  those  manufactured  by  Kullite  Semiconductors 
Inc.  which  are  made  of  p-type  silicon.  These  gauges  are  evaporated 
on  glass  and  alumina  substrates  and  they  might  be  used  to  measure 
shear  waves  in  these  materials.  However,  their  shear  piezoresis¬ 
tivity  will  probably  be  influenced  by  both  the  pressure  and 
temperature  rise  which  result  from  the  leading  compressional  wave. 
Thus,  more  calibration  work  has  to  be  done  on  these  gauges  before 
they  can  be  used  for  quantitative  measurements  in  shock  wave 
experiments . 
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The  problem  of  gauge  slippage  because  of  weak  epO):y 
bonding  between  two  solid  specimen  disks  has  to  be  carefully 
analyzed  and  one  should  use  the  various  techniques  which  have  been 
developed  for  the  prevention  of  such  slippage.  The  two  mo3t 
promising  one  are  (1)  the  use  of  alumina-filled  epoxy  and  (2) 
embedding  the  gauge  in  precisely  prepared  grooves  in  the  target, 
without  using  epoxy  and  relying  on  the  friction  between  the  disks. 
The  problem  of  the  epoxy  failure  to  transmit  high  shear  stresses  is 
absent  if  one  is  using  liquid  specimens.  Thus,  it  could  be  very 
advantageous  to  perform  oblique  impact  experiments  on  liquids  and 
measure  the  shear  stresses  in  these  materials. 

In  order  to  check  the  feasibility  of  the  proposed 
technique,  an  obi ique- Impact  experiment  wa3  performed  with  the  50 
mm  gas  gun  of  the  University  of  Dayton  Research  Institute.  This 
gun  does  not  have  a  keyway  arrangement  for  oblique  impact 
experiments.  Thus,  in  order  to  achieve  a  reasonable  planarity  of 
the  impact,  a  f ew  shots  were  conducted  in  which  it  was  determined 
that  the  amount  of  rotation  the  projectile  experiences  during 
launch  was  about  3°*  Henceforth,  we  place  the  projectile  in  a 
position  rotated  by  the  same  amount.  A  check  for  this  compensation 
is  obtained  by  observing  the  rise  time  of  the  compressional  wave  as 
recorded  by  the  gauge.  Although  the  compensation  was  not  perfect, 
the  rise  times  were  good  enough  for  a  qualitative  evaluation  of  the 
techni que . 

Next  a  commercial  mangaritn  gauge  was  placed  between 
two  98  percent  pure  alumina  disks.  This  material  had  a  density  of 
3-75  g/ cc  and  a  longitudinal  sound  velocity  of  =  10. *1  mm/ys. 
Based  on  past  work  the  target  response  should  be  elastic.  The 
thickness  of  the  first  disk  was  13.7  mm  so  that  the  compression  and 
shear  waves  should  have  a  large  separation  between  their  arrivals 
at  the  gauge  location.  The  impactor  was  a  thick  copper  disk,  wedge 
shaped  as  shown  in  Figure  25.  and  the  angle  of  obliquity  was  15°. 
The  impact  velocity  in  this  experiment  was  210  m/ s .  Using  these 
values,  together  with  an  estimated  shear  wave  velocity  of  b  mm/ys, 
the  shear  generated  at  the  impact  face  should  be  t»5  kbar. 

According  to  Reference  31,  this  shear  stress  is  larger  than  the 


maximum  shear  that  a  thin  layer  of  epoxy  can  transmit  dynamically 
(»3kb).  However,  the  intent  was  not  to  measure  the  shear  stress 
amplitude,  but  rather  to  record  its  arrival  at  gauge  location. 

The  gauge  record  in  this  experiment  is  shown  in 
Figure  26.  It  can  be  clearly  seen  that  the  arrival  of  two  waves  at 
the  gauge  location  occurs  with  a  separation  of  about  1.2  ys.  The 
rise  time  of  the  first  wave  is  relatively  large,  a  constant  stress 
plateau  is  not  observed  and  the  peak  stress  is  less  than  the 
predicted  stress,  as  calculated  from  the  impact  parameters. 
Nevertheless,  the  detection  of  two  separate  waves  reaching  the 
gauge  plane  is  unambiguous.  Using  the  measured  value  of  »  1 0 . M 
rnm/ps  (with  ultrasonic  techniques)  and  the  time  interval  between 
the  waves  we  get  a  value  of  C  -  5.4^  for  the  shear  wave  velocity. 
This  value  is  within  the  range  of  values  for  various  alumina 
specimens,  as  given  by  Coors,  with  a  range  from  5.06  to  6.38  rara/us. 

The  relatively  large  apparent  shear  wave  amplitude 
does  not  necessarily  mean  that  the  coefficient  of  manganin  is 

large.  The  more  probable  explanation  is  that  the  epoxy  layer 
between  the  two  alumina  disks  could  not  transmit  the  large  shear 
stress  (of  5  kb)  and  it  slipped,  causing  the  gauge  to  elongate  and 
its  resistance  to  increase.  This  can  be  verified  by  performing 
more  experiments  with  lower  shear  stresses  and  with  gauges  oriented 
differently.  These  different  orientations  will  result  in  different 
elongations  when  slippage  occurs,  thus,  different  resistance 
changes . 

In  summary,  we  have  proposed  to  use  piezoresistance 
gauges  to  monitor  both  compression  and  shear  waves  in  oblique- 
impact  experiments.  The  first  step  is  to  identify  and  statically 
calibrate  suitable  gauge  materials  that  possess  large  shear 
piezoresistance  coefficients.  The  second  step  is  to  embed  the.se 
gauges  in  well-defined  elastic  solids,  like  ceramics,  in  order  to 
demonstrate  their  performance  and  check  their  static  calibration. 
After  these  two  steps,  they  can  he  embedded  in  any  specimen  using 
strong  epoxy  adhesive  (as  suggestea  in  Rex".  3i  )  or  the  fine 
grooving  technique  which  was  used  in  Ref.  8.  We  suggest  Ytterbium 
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Figure  26.  Manganin  Gage  Record  from  the  Oblique-Impact 
Experiment,  [Note  the  arrival  of  the  two 
waves.  (0.5  V/div;  lys/div.)] 
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gauges  because  this  material  has  large  and  ^  coefficients  and 

it  may  have  a  large  coefficient.  The  experiments  should  be 

performed  in  a  properly  keyed  gun  for  good  planarity  of  the  impact. 
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SECTION  III 
ANALYTICAL  MODELING 


Elementary  constitutive  models  relate  stress  to  strain  and 
are  generally  independent  of  loading  rate  and  load  history.  While 
these  are  easy  to  use,  they  are  not  accurate  for  many  metals  under 
load  cycles  of  engineering  importance.  In  particular,  numerical 
simulations  of  projectile  penetration  and  explosive  forming  require 
more  accurate  constitutive  models  than  are  usually  available. 
Rate-dependency  is  taken  into  account  by  overstress  concept  based 
viscoplastic  models,  for  example,  those  of  Perzyna  [32]  and  Malvern 
[33].  However,  these  models  do  not  consider  loading  history.  A 
state  variable  based  viscoplastic  model  must  be  employed  in  order 
to  describe  both  rate-  and  history-dependent  behavior. 

The  state  variable  model  of  Bodner  and  Partom  [34]  (B-P)  has 
many  attractive  features:  (1)  ability  to  model  a  broad  range  of 
load  histories,  (2)  ease  of  determination  of  parameter  values,  and 
(3)  adaptability  to  finite  difference  and  finite  element  computer 
codes.  This  model  has  been  successfully  used  by  various  inves¬ 
tigators  to  describe  several  structural  alloys  under  cyclic  and 
time  dependent  loading  conditions.  The  strain  rates  Involved  in 

— *  P  —  1 

those  studies  were  between  10  and  10s  .  The  present  work 

reports  the  successful  application  of  the  B-P  model  to  material 
behavior  at  high  strain  rates  (>100  s  ).  A  unique  algorithm  was 
developed  to  evaluate  the  model  parameters  from  experimental  data. 

Details  of  the  modeling  technique  have  been  presented  in  an 
earlier  report  [1],  For  completeness,  we  describe  briefly  the 
Bodner-Partom  model  and  the  model  parameter  evaluation  scheme.  We 
also  used  the  model  to  describe  high  temperature  SHB  test  results. 
The  following  section  discusses  these  various  results. 

3.1  BODNER-PARTOM  MODEL 

We  present  the  major  features  here,  as  well  as  an  extension 
of  the  model  for  extreme  work  hardening  materials,  as  proposed 
earlier  by  Bodner  et  al .  [35]. 
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The  total  strain  rate  is  assumed  to  be  decomposable  into 
elastic  and  Inelastic  components. 


e .  . 

ij 


•  e  •  p 

e .  .  +  et  . 
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(19) 


Both  are  nonzero  for  all  loading/unloading  conditions.  The  elastic 
strain  rate  is  related  to  the  stress  rate  by  the  elastic  constants 
(Hooke's  Law).  The  inelastic  strain  rate  is  assumed  to  be  a  func¬ 
tion  of  stress,  o^,  a  state  variable,  Z,  and  to  follow  the 
Prandtl-Reuss  flow  rule 


i  j 


1  J 


XS; 


1  J 


(20) 


where  e^j  and  S^j  are  the  deviatoric  plastic  strain  rate  and 
stress,  respectively.  Squaring  Equation  (20)  gives 
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where  is  the  second  invariant  of  the  plastic  strain  rate  and 
is  the  second  invariant  of  the  stress  deviator. 
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where  Dq  is  the  limiting  value  of  the  plastic  strain  rate  in  shear, 
and  n  is  a  parameter  that  is  mainly  related  to  strain  rate  sen¬ 
sitivity.  Z  is  a  measure  of  the  overall  resistance  of  the  material 
to  plastic  flow,  and  it  depends  on  the  loading  history. 

The  plastic  strain  rate  can  be  expressed  in  the  following 
form  by  regrouping  Equations  (20)-(22) 
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It  is  assumed  that  Z 


i3  a  function  of  the  plastic  work, 
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1  -  m  (Z1  -  Z) 


(24) 


Integration  of  this  equation  will  lead  to: 

-mW 

Z  -  Z1  -  (Z1  -  Zo)e  P 


(25) 


where  W  ■=  o,  .ep  .  Z  is  the  Initial  value  of  Z,  Z1  is  the  maximum 
P  ij  ij  o  ' 

value  that  Z  can  attain,  and  m  is  a  parameter  that  describes  the 
strain  hardening  behavior  of  the  material. 

To  describe  extreme  strain  hardening  material  behavior,  m  can 
be  assumed  to  take  the  following  form: 


-aW 

m  =  m  +  m,  e  p  (26) 
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The  expression  for  Z  is  obtained  by  combining  Equations  (24)  and 
(26)  and  integrating, 
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The  additional  constants,  mQ  and  m^ ,  replaces  the  original  '  m' .  We 
also  need  'a'  to  model  extreme  strain  hardening  behavior  of 
material  like  annealed  OFHC  copper. 


The  effect  of  temperature  on  the  flow  stress  is  included 
through  the  parameter  n.  Based  on  experimental  evidence,  it  is 
assumed  to  vary  as  an  inverse  function  of  absolute  temperature  (T): 


n  =  A  +  y 


(28) 


where  A  and  B  are  model  parameters  and  T  is  expressed  in  degrees 
Kelvin.  ThU3  the  B-P  model  contains  five  principal  material 
parameters  DQ,  ZQ(  Z^  m,  and  n  that  have  to  be  evaluated  from  high 
strain-rate  experiments  under  room  temperature  for  a  material  with 
negligible  strain  hardening.  For  a  strain  hardening  material,  the 
number  of  constants  is  increased  by  two,  with  m  being  replaced  by 
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mQ ,  ,  and  a.  However,  to  extend  the  model  to  elevated  tempera¬ 

ture,  we  need  to  express  n  through  the  additional  parameters  A  and 
B. 


3.2  IMPLEMENTATION  INTO  STEALTH  CODE 


In  most  general  purpose  finite  difference  codes,  the  plastic 
flow  is  described  by  rate  independent,  Von  Mises-yield-criterion- 
baaed,  models.  For  instance,  in  the  STEALTH  [36]  finite  difference 
code  that  we  used  in  our  study,  the  stress  deviators  are  esti¬ 
mated  based  on  the  elastic  material  behavior  at  the  end  of  each 
time  step.  The  effective  stress,  oeff»  is  estimated  using: 
oeff  -  /3/S  STjS~j .  °eff  is  then  compared  with  the  constant  yield 
stress,  Yq  which  is  an  input  into  the  code.  When  oef{,  exceeds  YQ, 
the  stress  deviators  are  brought  back  to  the  yield  surface  by  a 
scaling  factor. 


The  current  algorithm  to  describe  the  B-P  model  requires  an 
iterative  scheme.  The  plastic  strain  rates  in  Equation  (23)  are 
calculated  initially  based  on  the  elastic  stress  deviators.  Since 
the  total  strain-rates  are  known  at  the  beginning  of  each  global 
time  step  of  the  finite  difference  solution  scheme,  we  can  compute 
based  on  the  initial  estimate  of  Using  Hooke’s  law,  the 

stress  rates  o.  .  are  estimated.  We  can  now  estimate  S,  .  based  on 
the  With  the  improved  values  of  3^,  the  plastic  strain  rates 

are  again  calculated.  This  procedure  is  continued  until  the  values 
of  5 ^ j  computed  for  the  two  successive  iterations  converged  to  the 
same  values  within  the  tolerant  limits.  We  introduced  a  sub- 
incremental  time  stepping  scheme  during  the  iterations.  The  global 
time  is  further  divided  into  small  steps.  The  solutions  converged 
within  2  or  3  iterations.  Rajendran  and  Grove  [37]  have  discussed 
this  technique  in  detail. 


3.3  ESTIMATION  OF  B-P  PARAMETERS 

There  are  seven  B-P  parameters,  D, .  n,  m_ ,  m, ,  a,  and  Z, . 
Rajendran  et  al .  described  the  estimation  scheme  in  detail  in 
Reference  38.  For  completeness,  a  brief  summary  of  the  scheme  that 
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was  employed  in  determining  these  parameters  is  given  in  this 

section.  D  in  the  model  is  the  maximum  strain  rate  which  the 
o  8  "  1 

material  can  experience;  it  was  arbitrarily  set  at  10s  .  The 

parameters  n,  ZQ,  Z1  and  mQ  were  first  estimated  from  tensile  SHB 

data,  and  then  modified  to  be  consistent  with  both  the  SHB  test  and 

the  plate  impact  data.  The  parameter  m1  was  initially  assumed  to 

be  zero,  and  only  in  the  case  of  OFHC  copper  [1,26]  does  it  have  a 

nonzero  value. 

A  series  of  interactive  automated  computer  programs  was 
developed  to  carry  out  the  evaluation  procedure.  A  program  called 
BPSOLVE  [39]  was  written  to  generate  uniaxial  stress  loading 
curves.  The  program  was  used  to  estimate  n,  mQ,  Zq  and  Z 1  from  the 
SHB  tension  test  results.  Next,  STEALTH  was  used  to  simulate  plate 
impact  tests;  usually  the  value  of  n  was  adjusted  to  match  the 
yield  behavior.  The  SHB  tests  were  resimulated  with  the  new  value 
of  n,  and  the  other  parameters  were  adjusted  interactively.  If 
necessary,  further  refinements  in  n  were  then  made  to  maintain 
adequate  agreement  with  the  plate  impact  data. 

Since  the  temperature  dependence  is  modeled  through  the 
parameter  n,  the  other-  parameter  values  were  assumed  to  be  the  same 
under  elevated  temperatures.  Several  SHB  tests  were  conducted  at 
different  temperatures.  Using  BPSOLVE,  the  value  of  n  correspond¬ 
ing  to  each  temperature  was  determined.  We  then  constructed  an  n 
vs.  1 /T  plot.  From  this  plot,  A  and  B  were  determined. 

We  also  present  the  new  Bodner  model  [40]  which  ha3  been  used 
in  our  study  to  describe  OFHC  copper.  The  following  section  very 
briefly  discusses  the  one-dimensional  version  of  this  model. 

3.4  BODNER  MODEL  FOR  DIRECTIONAL  HARDENING 

In  this  model,  the  state  variable  Z  is  assumed  to  be  made  up 
of  both  isotropic  and  directional  components, 

Z  -  Z1  +  ZD  (29) 

The  evolution  equation  for  Z1  is  given  by: 
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(30) 


21  -  m1  (Z,  -  z1)  wp 

with  the  initial  condition,  Z^O)  -  Z  ,  Z1  is  the  limiting 

I  0 

(saturation)  value  of  Z  ,  m.  is  the  hardening  rate.  The  expression 
D 

for  Z  is  given  by: 


Z 


D 


(31  ) 


where  u^  are  the  direction  cosines  of  the  current  stress  state. 
Under  uniaxial  stress  state, 


"  B  |~g'|~  where  |o| 
evolution  equation  for  B  is, 


is  the  absolute  value  of  o. 
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This  new  Bodner  model  assumes  a  simplistic  plastic  strain-rate 
equation  as: 
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The  term  (n+1/n)  that  appeared  in  Equation  (22)  of  the  B-P  model  is 
not  used  in  the  new  formulation.  This  model  has  the  same  number  of 
parameters  as  the  B-P  model  for  high  strain  hardening  materials,  as 
in  the  case  of  OFHC  copper.  The  parameters  are:  DQ,  n,  m1 ,  m2,  ZQ , 
Zy  and  Z ^ .  We  used  this  model  to  describe  the  annealed  OFHC 
copper.  The  SHB  test  results  reported  in  Reference  1  were  used  in 
the  modeling  efforts. 


3.5  MODEL  PARAMETERS 

We  reported  the  B-P  model  parameters  for  1020,  C1008,  HY100 
steels,  OFHC  copper  and  7039-T611  aluminum  in  an  earlier  report  Cl  3. 
In  "Edition  to  these  five  materials,  we  determined  model  parameters 
f.,.  6061  -T6  aluminum,  Nickel  200,  W-2  Tungsten,  and  Armco  Iron.  We 

also  extended  the  model  parameters  for  C1008  steel  to  include 
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TABLE  6 


BODNER-PARTOM  PARAMETERS 


MATERIAL 

(kbar) 

(kblr) 

- If - 

n 

(kbar 

A 

) 

B 

Cl 008  Steel 

55 

70 

0.  4 

1.5 

0.245 

46 

HY100  Steel 

24 

35.5 

1  .  2 

1 

- 

- 

1020  Steel 

6.4 

9.3 

4.0 

3.0 

- 

- 

6061  -T6 

Aluminum 

4.5 

5.5 

4.0 

12.0 

-2.86 

2343.24 

6039-T64 

Aluminum 

5.6 

7.6 

4.0 

2.8 

- 

Nickel  200 

3.2 

8.  2 

4.0 

4.0 

- 

- 

W2-  Tungsten 

87.5 

1  00.0 

0.578 

15.0 

0.1  66 

134.3 

Armco  Iron 

26.46 

42.00 

0.58 

5.65 

- 

- 

T 

OFHC  Copper 
(annealed ) 

8.0 

65.5 

0.4 

1  .  4 

*  For  high  temperatures,  value  for 
n  =  A  ♦  B/T. 

n  is  calculated 

based  on 

^  Two  additional 
a»1  50  kbar 

constants 

for  OFHC 

copper 

are: 

=  15  kbar  ^ 

and 
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temperature  effects.  An  attempt  was  al30  made  to  improve  the 
modeling  of  OFHC  copper  through  the  new  Bodner  model  which  con¬ 
siders  directional  hardening  effects.  Table  6  provides  the  B-P 
model  parameters  for  the  various  metals.  The  model  simulations  of 
SHB  and  plate  impact  results  are  shown  in  Figures  27-35.  It  can  be 
seen  from  these  figures  that  the  model  reproduced  the  experimental 
data  extremely  well. 

3.5.1  Improved  Model  for  OFHC  Copper 

Figure  36  shows  the  SHB  test  data  at  e  =  1100s  ,  and 

the  corresponding  simulation  of  B-P  model.  The  model  simulated  the 
test  results  reasonably  well.  However,  the  unrealistic  shape 
(between  points  A  and  B  in  Figure  36)  of  the  simulated  curve  indi¬ 
cates  the  inadequacy  of  B-P  mooel  to  handle  large  strain  hardening 
behavior . 

Analysis  (see  Reference  1)  of  the  B-P  model  equations 
showed  that  this  inadequacy  cannot  be  eliminated  due  to  the  B-P 
model  formulation  baaed  on  plastic  work.  V,'c  attempted  to  model 
OFHC  copper  using  the  new  Bodner  model.  The  model  parameters  were 
obtained  using  the  BPSOLVE  program.  It  required  a  slight  modifica¬ 
tion  of  BPSOLVE.  The  corresponding  constants  are  given  in  Table  7. 

TABLE  7 

BODNER  MODEL  PARAMETERS 


n 

^  -1 
(kbar  ) 

(kbar  1  ) 

Zo 

(kbar ) 

Z1 

( kbar ) 

Z3 

( kbar ) 

0.3*1 

2.1 

1  50 

1  0 

370 

1  1  0 

Figure  37,  shows  the  actual  SHB  test  results  with  the 
B-P  model  and  Bodner  model  simulations.  It  can  be  seen  from  the 
Figures  38-39  that  the  concavities  in  the  simulated  stress-strain 
curves  at  three  different  strain  rates  could  not  be  eliminated  even 
with  the  new  model.  The  model  simulations  are  also  separately 
plotted  in  these  figures.  The  concavity  is  not  that  much 


Figure  27.  Comparisons  of  B-P  Model  Simulations  with  SHB 
Room  and  High  Temperature  Tests  Results  for 
C1008  Steel. 
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Figure  27.  Comparisons  of  B--P  Model  Simulations  with  SHB 
Room  and  High  Ternoerature  Tests  Results  for 
C100S  Steel  (concluded). 
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Figure  28.  Comparison  of  B-P  Model  Simulation  with  Plate 
Impact  VISAR  Data  for  C1008  Steel. 
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Figure  29.  Comparison  of  B-P  Model  Simulations  with  SHB  Room  and  High  Temperature  Tests 
Results  for  6061-T6  Aluminum. 
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Figure  30.  Comparison  of  Quasi-static  Tensile  Test  Results  with  B-P  Model  Simulations 
for  W-2  Tungsten. 
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Figure  31.  Comparison  of  B-P  Model  Simulations  with  SHB  Room  and  High  Temperature 
Test  Results  for  W-2  Tungsten  (concluded) . 
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Figure  32.  Comparison  of  B-P  Model  Simulation  with 
Manganin  Gauge  Stress  History  in  PMMA 
W- 2  Tungsten . 
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Figure  33.  Comparison  of  B~P  Model  Simulations  with 
SHB  Tensile  Test  Data  for  Nickel. 
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Figure  34.  Comparison  of  B-P  Model  Simulations  with 
SHB  Tensile  Test  Data  for  Armco  Iron 
(concluded) . 
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Comparison  of  B-P  Model  Simulations  with 
Plate  Impact  Test  Plot  for  Armco  Iron. 
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Figure  36.  Comparison  of  B-P  Model  Simulation  with  SHB 
Test  Data  for  OFHC  Copper.  (Note  the 
unrealistic  shape  between  points  A  and  B.) 
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Figure  38.  Comparisons  Between  B-P  Model,  New  Bodner  Model,  and  SHB  Tests  Data  for 
OFHC  Copper  at  e  =  1100s-1. 
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Figure  39.  Comparisons  Between  B-P  Model,  New  Bodner  Model,  and  SHB  Test  Data  for 
OFHC  Copper  at  £  =  1300s“l. 


pronounced  in  the  new  model  when  compared  to  the  B-P  model.  It 
further  confirms  that  the  plastic  work  based  evolution  equations 
for  state  variable  may  not  be  able  to  describe  extreme  strain 
hardening  behavior. 

As  a  final  exercise,  we  slightly  modified  the  evolu¬ 
tion  equation  for  Z  in  the  B-P  model  formulation  as  follows: 

2  -  m  (Z,  -  Z)  ePff  (34) 

•  n 

where  is  the  effective  plastic  strain  rate  and  is  defined  by, 


•  p 

f  r 

ef  f 


7 


"5  *p  rp 
3  ei jei j 


Under  uniaxial  stress  condition,  it  reduces  to: 


(35) 


eP  =  ep 

eef  f  1 


(36) 


We  modeled  the  ShB  test  results  on  UFHC  copper  using 
this  modified  equation.  Figure  40  shows  the  comparison  of  the 
model  simulation  with  the  SHB  stress- strai n  curve.  Use  of  effec¬ 
tive  plastic  strain  rate  has  eliminated  the  unrealistic  concavity 

in  the  stress- at rai n  curve.  It  also  compared  with  the  test  results 

•  •  n 

very  well.  In  Figure  41  the  model  simulations  of  Wp  and  based 

B-P  models  are  compared.  It  can  be  seen  from  this  figure,  the 

•  n 

shape  of  the  stress-strain  curve  of  based  model  is  extremely 

VI  l 

good.  The  strain-hardening  behavior  is  very  realistic  due  to  the 
absence  of  the  eoncave-up  behavior.  At  present  we  are  not  sure 
whether  this  indicates  that  the  B-P  model  equation  requires 
modification.  The  stress- strai n  behavior  based  on  the  strain  rate 
jump  test  results  on  OFHC  copper  was  modeled  by  the  Bodner-Partom 
equations  [41]  extremely  well.  The  merit  of  B-P  model  is  in  its 
ability  to  model  the  loading  history  effect  on  flow  stress.  It  is 
yet  to  be  seen  whether  a  similar  match  can  be  achieved  with  the 

^P  T,  n 

t-  ^  uaacu  u  i  itivuvi.  • 
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Comparison  of  Modified  B-P  Model  (Effective 
Plastic  Strain  Rate  Based)  with  SHB  Test 
Data  for  OFHC  Copper. 


Figure  40. 
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Figure  41.  Comparison  Between  Effective  Plastic  Strain 
Rate  Based  Model  and  Plastic  Work  Based 
B-P  Mode  1 . 
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3-6  PRECURSOR  DECAY  ANALYSIS 


The  phenomenon  of  elastic  precursor  decay  in  planar  impact 
experiments  has  been  a  subject  of  great  interest  to  researchers  in 
dynamic  plasticity  and  shock  wave  propagation  during  the  last  few 
decades.  Efforts  to  understand  and  predict  precursor  decay  have 
involved  both  experiments  and  analy t i cal/ compu t at i onal  modeling. 

The  latter  has  been  largely  associated  with  concepts  based  on 
dislocation  dynamics.  In  this  investigation,  we  invoke  the  basic 
assumption  that  the  precursor  decay  phenomenon  is  associated  with 
rate-dependent  material  response.  Our  objective  was  to  model  this 
phenomenon  using  the  Bodner-Partom  model.  In  the  present  study, 
the  B-P  model  parameters  for  rate-dependent  C1008  steel  were  deter¬ 
mined  based  on  a  combination  of  dynamic  stress-strain  curves 
obtained  from  the  SHB  tests  and  the  HEL  obtained  from  planar  plate 
impact  tests.  Using  the  STEALTH  code,  we  simulated  plate  impact 
tests.  The  transient  wave  profiles  in  terms  of  stress  vs.  time  and 
velocity  vs.  time  were  numerically  obtained.  The  B-P  model  predic¬ 
tions  of  these  wave  profiles  were  compared  directly  with  the 
experimental  data  reported  in  Reference  1.  In  the  numerical 
evaluation,  a  mathematical  definition  of  the  HEL  has  also  been 
proposed . 

Precursor  decay  is  a  phenomenon  that  is  associated  with  a 
planar  plate  impact  test  configuration.  It  is  defined  as  the  decay 
of  the  longitudinal  3tress  amplitude  (o^^)  of  the  elastic  wave 
(precursor)  from  a  peak  value  at  the  plane  of  impact  (-  1/2pcv,  p, 
the  density,  c,  the  sound  speed,  v,  the  impact  velocity)  to  a 
steady  state  value  at  some  distance  from  the  impact  plane.  For  an 
elastic-perf ectly  plastic,  rate  independent  material,  the  oHF[  is 
related  to  the  initial  yield  stress  (Y  )  under  uniaxial  stress 
through  the  relationship  (5),  given  in  Section  2.^. 

Precursor  decay  cannot  be  modeled  by  rate  independent  models. 
A  rate  dependent  theory  is  required  to  model  precursor  decay.  When 
such  a  theory  is  employed  in  a  finite  difference  code  to  evaluate 
the  decay  phenomenon,  a  mathematical  definition  o  1'  HEL  is 
essential.  For  this  purpose,  using  a  concept  similar  to  the  0.2 


91 


percent  offset  criterion  that  has  been  used  to  determine  the  in¬ 
itial  yield  stress  in  a  uniaxial  quasi-static  tension  test,  a 
critical  value  of  effective  plastic  strain,  is  proposed.  A 

detailed  analysis  of  the  proposed  criterion  can  be  found  in 
Ref erence  42 . 

To  determine  the  critical  value  of  ^ ,  we  simulated  a  plate 
impact  experiment  in  which  the  flyer  and  target  plates  were  3  mm 
and  20  mm,  respectively.  The  precursor  decay  was  calculated  first 
using  a  critical  value  of  0.002  and  then  a  value  of  0.0002.  Figure 
42  compares  the  decay  results  obtained  in  these  two  cases.  The 
steady  state  HEL  value  for  C1008  steel  obtained  from  plate  impact 
experiments  [1]  was  around  9  Kbar  and  this  compared  well  for  - 

0.0002  as  can  be  3een  from  Figure  42.  We  also  employed  a  value  of 
0.0001  which  did  not  significantly  change  the  results.  The  value 

of  0.0002  was  thus  arbitrarily  chosen  as  the  critical  value  of 
eP 
ef  f  ‘ 

The  experimental  HEL  points  were  compared  in  Figure  43  with 
the  results  obtained  from  STEALTH  using  the  B-P  model  parameters. 
The  analytical  prediction  of  precursor  decay  matched  extremely  well 
with  the  experimental  results.  The  proposed  definition  of  HEL  was 
also  verified  by  employing  the  viscoplastic  model  of  Perzyna  L  32  3  - 
We  validated  further  the  proposed  mathematical  definition  of  HEL  by 
modeling  precursor  decay  in  aluminum  7039-T64.  In  all  cases,  a 
critical  value  of  0.0002  for  ej^j,  was  employed  successfully. 

The  present  study  demonstrated  that  (a)  the  precursor  decay 
phenomenon  can  be  predicted  using  a  state  variable  based  viscoplas¬ 
tic  constitutive  model,  (b)  the  model  parameters  need  not  be 
calibrated  using  the  precursor  decay  data,  and  (c)  a  mathematical 
definition  of  HEL  is  required  to  model  the  decay  phenomenon.  We 
also  showed  the  ability  of  the  Bodner-Partom  model  in  modeling 
material  response  under  a  wide  range  of  stress  states  and  strain 
rates.  The  model  parameters  that  were  determined  using  SHB  and  the 
steady  state  value  of  HEL  are  general . 
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Figure  42. 


Precursor  Decay  along  the  Target  Thickness 
Using  Criterion. 
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COMPARISON  WITH  EXPERIMENTS 


Figure  43,  Comparison  of  the  Experimental  Data  with  the 
Model  Prediction. 
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SECTION  IV 

DYNAMIC  FAILURE  MODEL 


Dynamic  failure  processes  in  ductile  metals  has  been  a  study 
of  interest  for  the  last  few  decades.  There  are  numerous  journal 
articles,  technical  reports,  and  conference  proceedings  on  this 
subject.  Investigators  with  widely  varied  background,  such  as 
applied  mechanists,  metallurgists,  shock  physicists,  etc,  have  been 
contributing  to  the  understanding  of  dynamic  fracture.  The  fun¬ 
damental  failure  mechanism  associated  with  dynamic  fracture  in 
ductile  metals  is  one  which  considers  the  failure  process  as  being 
initialed  by  the  nucleation  of  voids  around  inclusions  and  their 
subsequent  growth  and  coalescence  as  suggested  by  MeClintock  [43]. 
During  the  last  two  decades,  different  approaches  have  been  taken 
by  several  investigators  to  apply  this  mechanism  as  a  means  of 
predicting  ductile  failure.  Among  them,  Rice  and  Tracey  [44] 
considered  the  growth  kinetics  of  a  single  void  in  the  matrix  and 
developed  a  model  which  showed  the  dependence  of  void  growth  rate 
on  the  triaxiality  of  the  stress  state. 

Hancock  and  Mackenzie  [4b]  provided  the  experimental  evidence 
to  support  this  dependency  on  the  stress  state  and  used  this  model 
to  predict  failure  initiation  in  notched  tensile  specimens  under 
quasi-static  loading  conditions.  Their  study  is  under  obviously 
low  mean  stress  levels. 

There  are  other  experimental  configurations  in  which  this 
fundamental  failure  mechanism  is  frequently  observed  and  reported. 
One  such  experiment  involves  impacting  a  flyer  plate  of  one 
material  against  a  target  of  the  same  or  other  material  at  a  higher 
velocity.  Compressive  stresses  are  produced  and  transmitted  im¬ 
mediately  from  the  plane  of  impact  to  the  adjacent  stress-free 
areas  of  the  material  in  the  form  of  a  stress  pulse.  Fracture  is 
then  induced  by  tension  arising  from  the  interaction  of  reflected 
shock  waves  from  the  stress  free-planes  parallel  to  the  plane  of 
impact.  The  interesting  feature  is  the  growth  and  coalescence  of 
voids  under  high  mean  stress  levels  leading  to  failure  in  the 
target  plate;  this  type  of  failure  is  often  referred  to  as 
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spallation.  Earlier  spall  models  were  mostly  based  on  a  simple 
criterion  in  which  stress,  strain,  pressure,  or  energy  is  assumed 
to  reach  a  critical  value.  Failure  had  been  assumed  to  occur 
instantaneously.  This  type  or  time  independent  failure  model  may 
work  well  when  the  spall  is  well  above  the  threshold  conditions. 
Later,  Tuler  and  Butcher  [46]  introduced  a  time  dependent  failure 
criterion.  They  assumed  that  failure  occurs  only  after  a  critical 
value  of  damage  is  reached.  This  model  has  three  parameters  which 
have  to  be  determined  from  spall  experiments.  Several  other  inves¬ 
tigators  [47-50]  proposed  models  to  predict  the  spall  threshold. 
They  considered  the  nucleation  sad  growth  rates  of  void3  through 
rigorous  mathematical  modeling  of  the  failure  process.  Meyers  and 
Aimone  [51]  have  reviewed  the  fundamental  aspects  of  spallation  and 
the  modeling  efforts  by  the  various  research  workers. 

The  various  approaches  can  be  grouped  into  three  categories: 
(1)  micro-statistical  approach,  (2)  physically  motivated  empirical 
approach,  and  (3)  continuum  mechanics  based  approach.  The  first 
approach  was  mainly  pursued  by  Curran  and  co-workers  [52],  Based 
on  the  extensive  micro-statistical  analyses  of  the  void  sizes  and 
orientations,  they  developed  models  to  describe  the  void  nucleation 
and  growth  processes.  In  the  second  ad  hoc  approach,  the  objective 
was  to  develop  models,  not  based  on  any  extensive  analyses  of 
microstructures  or  the  kinetics  of  the  failure  processes,  but 
mainly  based  on  the  ability  to  predict  certain  measured  quantities 
such  as  the  free  surface  velocity  of  the  target  plate.  Under  thi3 
approach  most  of  the  works  had  been  strictly  developed  and  limited 
to  rectilinear  motions.  In  most  cases,  the  generalization  of  such 
theories  to  any  other  geometry  is  not  possible.  The  third  approach 
considered  a  general  solution  to  the  problem  of  spallation. 

Davison,  Stevens,  and  Kipp  [49]  presented  a  theory  of  spall  damage 
based  on  a  unified  and  thermodynamically  consistent  treatment  of 
e 1  as t i c- v i scoplas ti c  deformation.  A  continuum  mechanics  approach 
was  employed  with  a  three-dimensional  theory. 
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MODEL  FORMULATION 


1 

There  are  at.  least  four  phases  in  the  constitutive  modeling 
of  damaged  materials  that  were  initially  intact.  The  first  phase 
of  the  modeling  is  the  description  of  the  intact  material  behavior. 
Prior  to  the  development  of  damage,  the  aggregate  is  the  fully 
dense  matrix  material  which  is  usually  modeled  by  incompressible 
plastic  flow  theories.  The  second  phase  is  the  description  of  the 
damaged  material.  We  require  a  model  to  describe  the  behavior  of 
the  aggregate  material  wmch  contains  microvoids/microcracks.  In 
the  third  phase,  depending  on  the  nucleation  mechanism,  a  mathe¬ 
matical  description  of  the  process  will  be  required.  This  phase 
will  also  require  a  model  to  describe  the  growth  of  damage.  The 
last  phase  of  modeling  is  the  coalescence  of  damage  leading  to 
total  failure.  The  final  process  occurs  quite  rapidly  when  com¬ 
pared  to  the  nucleation  and  growth  processes.  Modeling  of  the 
fundamentally  observed  instabilities  [53.54]  during  the  final 
failure  process  will  be  cumbersome  and  tedious.  A  simplistic 
approach  is  to  describe  this  process  through  a  critical  criterion 
based  on  a  single  parameter.  For  example,  one  such  criterion  could 
be  based  on  a  critical  level  of  damage.  This  section  describes  the 
modeling  efforts  on  each  of  the  above  mentioned  phases. 

4,1.1  Phase  I  -  Modeling  of  Matrix  Material 

The  fully  dense,  void-free  matrix  material  can  be 
modeled  through  the  state  variable  baaeu  viseoplastic  constitutive 
equations  of  Boaner-Partom  [34].  The  B-P  model  will  be  used  in  the 
constitutive  model  formulation  for  the  porous  aggregate  in  terms  of 
equivalent  plastic  strain  rate  Dp  and  effective  stress,  Y^.  The 
equivalent  plastic  strain  rate  is  defined  by, 


Dp  -  VI  eP  ,ep  and  Y  -  /3J, 
m  *  3  i  J  i  J  m  2 


(37) 


i he  equation  for  Dp  can  be  obtained  from  equation  (??)  as 

m 


_  ,  ..  2n 

2  ~  _ i  inM  i  /  L 


Dm  ‘  -%  Do  exP'-l-2H>  1 

✓  3  m 


(38) 
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B-P  model  parameters  for  several  metals  have  been  already  deter¬ 
mined  and  reported  in  Section  III.  The  number  of  parameters  for 
describing  matrix  material  behavior  is  only  four,  i.e.,  m,  n,  2Q , 
and  Z  j  . 

Until  voids  nucleate,  the  aggregate  behavior  can  be 
described  by  the  B-P  model.  The  plastic  flow  in  the  void-free 
aggregate  is  incompressible,  i.e.,  the  sum  of  the  principal  or 
orthoganal  plastic  strains  is  equal  to  zero.  However,  the  nuclea- 
tion  of  voids  will  introduce  dilatation  and  the  plastic  yield 
behavior  will  depend  not  only  on  the  second  invariant  ,  but  also 
on  the  mean  stress  or  pressure.  The  constitutive  model  for  the 
aggregate  must  include  these  effects.  For  this  purpose,  we 
selected  a  yield-criterion-based  plastic  flow  rule  in  which  the 
pressure  dependence  enters  explicitly  into  the  calculations.  The 
following  section  describes  the  Phase  II  of  the  constitutive  model 
formulation  for  the  porous  aggregate. 

4.1.2  Phase  II  -  Modeling  of  Aggregate  with  Voids 

We  considered  a  yield-criterion-based  approach  in  the 
constitutive  model  formulation.  For  randomly  distributed  voids  or 
microcracks  contained  in  the  aggregate,  the  yield  behavior  will  be 
influenced  by  not  only  the  second  invariant  of  the  deviatoric 
stress  (J2)  but  also  by  the  pressure  or  mean  stress  (1^).  The 
following  form  of  the  yield  function  has  been  considered: 

A  (  p )  J2  +  B  (  p  )  I*  -  6  (  p  )  (39) 

where  A,  B,  and  <5  are  functions  of  relative  density,  p.  Ym  is  the 
effective  stress  in  the  material.  (Note:  the  subscript  ’ m'  means 
matrix  material  and  not  a  tensorial  index.)  Based  on  a  critical 
total  deformation  energy,  Doraivelu  et  al .  [55]  derived  the  follow¬ 
ing  expressions  for  A  and  B: 

A  (  p )  -  2  +  pP  (40) 

and 
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B (  p) 


1 


I 


( 4l  ) 


2 


The  expression  for  6 ( p )  was  given  by 

«(p)  -  2p2  -  1 


(42) 


This  particular  form  for  6 ( p )  was  chosen  based  on  the 
experimental  results  obtained  from  sintered  or  hot  pressed  powder- 
metallurgy  materials  under  uniaxial  compression  at  quasi-static 
strain  rates.  In  general,  this  function  is  material  dependent 
while  the  functions  A  and  B  are  independent  of  the  matrix  material 
behavior.  Thus,  the  yield  condition  for  the  aggregate  can  be 
written  as: 


*  -  (2  +  p2)  J2  +  I,  ~  6  (  p)  Y2  «  0 


(43) 


The  viscoplastic  strain  rates  in  the  aggregate  can  be 
calculated  using  the  flow  rule  derived  as: 


(44) 


The  proportionality  factor  A  can  be  obtained  using  the  flow  rule  in 
conjunction  with  the  following  relationship: 
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where  f  is  the  void  volume  fraction  and  related  to  relative  density 
'  p'  through  f  -  1-p.  The  above  expression  was  derived  from  the 
definition  that  the  total  plastic  work  in  the  aggregate  is  entirely 
due  to  the  plastic  work  done  by  the  matrix.  By  combining  equations 
(44)  and  (45),  X  can  be  expressed  by, 
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The  plastic  strains  in  the  aggregate  ean  be  written  as: 
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In  the  constitutive  model  formulation,  the  total  strain  rate  is 

decomposed  into  elastic  and  plastic  strain  rates.  The  elastic 

strain  rates,  e,  .  are  related  to  the  stress  through  Hooke's  law  as 
i  j 
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•  ,  o,  ■ 

ik  kj 


(48) 


where  D.,  is  the  inverse  of  elastic  modulus  matrix,  E, .  . 

l  k  ’ik 

Using  the  consistency  condition  which  holds  during 

the  plastic  flow,  we  can  obtain  an  expression  for  as: 
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An  expression  for  can  be  obtained  from  Equation  (48)  by  replac 

ing  the  elastic  strain  rate  as  the  difference  between  total  and 
plastic  strain  rates, 
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We  now  need  an  evolution  equation  to  calculate  f  at 
any  given  instance  of  the  loading  history.  The  void  volume  frac¬ 
tion  rate,  f  contains  two  parts;  one  due  to  the  growth  of  existing 
voids,  and  the  other  due  to  the  nucleation  of  new  voids. 

4.1.3  Phase  III  -  Nucleation  and  Growth  Models 


The  growth  law  can  be  directly  related  to  the  dilata 
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•  n 

where  repeated  index  means  summation  and  are  plastic  strain 

rates  in  the  three  principal  directions. 


The  most  widely  used  void  nucleation  model  was  the 
one  that  was  initially  used  by  Needleman  and  Chu  [56]  in  their 
analysis  of  localized  necking  in  biaxially  stretched  sheets.  The 
model  was  based  on  a  mechanism  in  which  voids  are  nucleated  due  to 
debonding  of  inclusions  from  the  matrix.  The  debonding  can  occur 
due  to  either  a  stress  or  a  strain  criterion  or  both.  The  cor¬ 
responding  model  is  given  by, 
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If  the  nucleation  is  due  to  only  the  matrix  debonding 
from  inclusions,  then  the  total  void  volume  nucleated  must  be 
consistent  with  the  volume  fraction  of  second  phase  particles. 
Therefore,  the  values  determined  for  the  parameters  f1  and  f2  must 
meet  this  requirement.  and  are  the  mean  equivalent  stress 

and  strain,  respectively  around  which  the  nucleation  stress  and 
strain  are  distributed  in  a  gausian  manner.  s^  and  s2  are  the 
standard  deviations  of  the  distributions.  These  two  parameters 
will  control  the  ranges  of  stress  or  strain  over  which  most  of  the 
voids  can  be  nucleated. 

The  evolution  equation  for  the  void  volume  fraction 
is  then  given  by  the  sum  of  f  and  f  as: 

O  1 1 
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where  Fq  and  F^  are  given  by  Equations  (53)  and  (5*0.  Six  model 
parameters  have  to  be  determined  to  describe  the  void  nucleation 
process,  when  the  process  is  due  to  both  stress  and  strain 
criteria. 

When  the  stress  state  is  triaxial  with  (P/Ym)  >>  1, 
the  nucleation  process  seems  to  be  stress-controlled;  whereas  under 
uniaxial  stress  states,  such  as  in  a  thin  ring  cylinder  or  under 
plane  stress  condition,  the  mechanism  seems  to  be  strain- 
controlled.  It  is  fairly  well  established  that  under  high  strain 
rate  and  high  pressure  loading  conditions,  the  void  nucleation 
mechanism  is  dominated  by  the  pressure  (mean  stress).  To  model  the 
spall  type  failure,  void  nucleation  due  to  F£  in  Equation  (55)  is 
negligible  and  in  turn,  we  can  set  Fe  »  0  in  the  calculations. 
Effectively,  the  void  volume  fraction  model  will  include  only  tnree 
parameters  which  can  be  determined  from  the  plate  impact  tests 
data . 

Equations  (38),  (47),  (99),  (50),  and  (55)  form  a  set 
of  relationships  that  has  to  be  solved  numerically  to  describe  the 
s tres3- strai n  response  of  a  void  contained  aggregate  material.  We 
developed  the  corresponding  numerical  algorithm  to  describe  these 
equations  in  the  STEALTH  finite  difference  code. 


9.2 


ANALYSES 


To  describe  spallation  under  a  one-dimensional  strain  impact 
condition,  we  modeled  a  plate  impact  test  #678  in  which  a  3  mm  1020 
steel  projectile  plate  impacts  a  6  mm  C1008  steel  target  at  a 
velocity  of  178  m/s.  We  did  not  consider  any  void  nucleation  due 
to  plastic  strain.  The  entire  nucleation  process  was  described  by 
the  stress  threshold  based  model.  The  nucleation  model  parameters 
that  were  to  be  determined  were:  s,  0^,  and  fQ.  Since  the  yield 
function  and  the  growth  law  did  not  contain  any  model  parameters, 
these  are  the  only  three  additional  parameters  that  need  to  be 
determined  to  describe  the  dynamic  failure  process. 
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Using  the  STEALTH  code,  we  modeled  the  plate  impact  test, 

#678.  The  finite  difference  simulations  of  this  test  with  various 
values  for  f  ,  and  s,  failed  to  generate  the  observed  spall 

signal.  The  calculated  velocity  vs.  time  (solid  line)  at  the  back 
free  surface  of  the  target  is  shown  in  Figure  t  *la .  The  experimen¬ 
tal  spall  signal  (dotted  line)  could  not  be  simulated.  It  can  be 
explained  by  analyzing  the  calculated  stress  history  in  the  spall 
plane  as  shown  in  Figure  Mb.  The  points  A-E  on  the  stress-time 
plot  correspond  to  the  points  A-E  on  the  veloci ty- time  plot.  The 
point  A  represents  the  Hugoniot  Elastic  Limit,  At  point  B, 

the  shock  wave  reaches  its  maximum  amplitude.  The  release  wave 
from  the  flyer  unloads  in  compression  at  point  C  and  reaches  its 
maximum  tensile  stress  at  point  D.  However,  in  the  stress-time 
plot,  the  tensile  amplitude  (10  kbars)  at  D  is  much  lower  than  the 
shock  (compression)  amplitude  (32  kbar)  due  to  stress  relaxation 
(due  to  voids  nucleation).  Note  that  a  value  of  20  kbars  is  as¬ 
sumed  for  oN.  Since  s  =  5,  the  nucleation  occurs  at  around  10-15 
kbars.  Due  to  void  growth,  the  stress  level  reduces  between  the 
points  D  and  E.  This  did  not  produce  any  spall.  Instead,  it 
created  a  porous  material  in  which  the  wave  speed  is  small  compared 
to  the  speed  in  the  intact  material.  This  is  the  reason  for  a 
longer  pulse  duration  under  the  tensile  loading.  In  the  velocity 
vs.  time  plot,  the  spall  nucleation  can  be  seen  at  point  D.  Since 
the  material  in  the  potential  spall  plane  did  not  loose  its 
strength  entirely,  the  spall  formation  is  not  complete.  It  is  also 
far  below  the  threshold  level  and  therefore,  the  velocity  continues 
to  decrease  between  points  D  and  E  in  Figure  Ma, 

If  we  introduce  a  spall  criterion,  such  as  o  -  ocr  at  point 
D,  a  spall  signal  could  have  been  created  as  shown  by  the  dotted 
line.  However,  our  objective  is  to  avoid  such  abrupt  failure 
criterion.  Instead  we  wanted  to  nucleate  and  grow  voids  until  the 
material  point  looses  its  strength  completely.  Theoretically,  a 
value  of  f  -  1  will  fail  the  material  completely. 

The  main  problem  that  we  faced  with  the  propose  d  i'  o  IT  m  \a  lat i on 
is  due  to  the  growth  law.  Since,  the  growth  model  did  not  contain 
any  model  parameter,  the  amount  of  void  growth  solely  depends  on 
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che  form  of  the  yield  criterion.  As  we  mentioned  eariier,  the  form 
of  <5  ( p )  in  the  yield  function  depends  on  the  material  and  loading 
regimes.  We  decided  to  modify  or  introduce  a  new  function  for 
<5(p).  Any  function  that  describes  6(p)  must  satisfy  two  basic 
conditions.  First,  at  p  »  1  (fully  dense  condition),  6(p)  =  1,  and 
secondly,  at  p  *  pcr,  6(p)  =  0,  where  pcr  is  a  critical  value  of  p 
at  which  the  material  looses  its  strength  completely.  Various 
forms  for  6 ( p )  were  considered  and  the  following  form  was  selected 
after  several  iterations, 


v  P 

6 ( p )  »  e  (56 ) 

where  B  is  a  material  parameter.  We  assumed,  pcr  -  0,  in  the  above 
function . 

We  introduced  this  equation  into  the  yield  function  repre¬ 
sented  by  equation  (43).  Figure  45  shows  the  effect  of  B  on  the 
yield  function  for  values  of  p  *  0.99,  0.85,  and  0.75.  In  terms  of 
void  volume  fraction  (f  °  1-p),  f  -  0.01,  0.15,  and  0.25.  The 

plots  of  /3 vs.  P/Y  show  that  the  amount  of  strength  loss  due 
2mm 

to  the  presence  of  voids  can  be  effectively  controlled  through  the 
parameter,  B. 

With  the  newly  introduced  function  for  6(p),  the  plate  impact 
test  on  C1008  steel  was  resimulated.  A  value  of  0.1  for  f  was 
assumed  in  our  simulations.  We  chose  a  standard  deviation  of  2  for 
s.  A  first  estimate  for  can  be  made  from  the  experimental 
velocity  vs.  time  plot.  We  can  assume, 


where  AVs  is  the  difference  between  the  measured  peak  particle 
velocity  and  the  particle  velocity  at  the  spall  signal.  This  is 
shown  in  Figure  46.  is  the  material  density,  and  C0  is  the 

sound  speed  (/(K+4/3G) 7p ) .  The  value  of  B  was  adjusted  until  the 
simulated  spall  signal  matched  well  with  the  experimental  data. 
Figure  48  compares  the  model  simulated  velocity  vs.  time  plot  with 
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Effect  of  the  Parameter  B  on  the  Yield  Surfaces. 


tne  experimental  result.  In  this  simulation,  the  following  values 
for  the  model  parameters  were  used. 

TABLE  8 

DYNAMIC  FAILURE  MODEL  PARAMETERS 


D 

0 

( 1  /sec ) 

m 

n  (kbar  ^) 

Z 

0 

(kbar) 

Z1 

(kbar ) 

3  f 

(kbar)  0 

°N 

(kbar ) 

B 

O 

CD 

0. 4  1.5 

55 

70 

2  0.  1 

21  .5 

20 

Since  the  spall  signal  of  test  #678  [1]  was  recorded  for  a 
short  time  duration  only,  the  model  simulation  was  compared  only  up 
to  point  G  as  shown  in  Figure  47.  The  slight  kink  at  point  F  was 
due  to  digitizing  of  the  experimental  data.  The  match  between  the 
points  A-E  was  due  to  the  Bodner-Partom  model  and  beyond  the  point 
E,  it  was  controlled  by  the  failure  model. 

To  further  validate  the  new  model,  a  plate  impact  test  (#038) 
in  which  a  copper  flyer  of  2  mm  thick  impacted  a  copper  target 
plate  of  9  mm  thick  at  a  velocity  of  185  m/s  was  simulated.  Using 
the  Bodner-Partom  model  constants  for  OFHC  copper  obtained  earlier 
by  Rajendran  and  Bless  [1]  and  given  in  Section  III  of  thi3  report, 
test  #538  was  simulated.  The  failure  model  parameters  that  could 
reproduce  the  experiment  are: 

3  °N 

(kbar)  (kbar) 

2  16.5 

The  velocity  vs,  time  plot  obtained  using  V1SAR  and  the  correspond¬ 
ing  model  simulation  plot  are  shown  in  Figure  48.  Since  OFHC 
copper  has  a  very  low  value  of  (<  1  kbar),  the  elastic  shock 

portion  of  the  velocity  plot  is  not  well  resolved.  Beyond  the  HEL, 
the  simulated  plot  compared  extremely  well  with  the  experimental 
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Figure  47.  Simulation  of  Test  #678  C1008  Steel  Using 

The  Failure  Model  with  New  Function  for  6(p). 
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Figure  48.  Comparison  of  Model  Simulation  with  the 
Experimental  Velocity  vs.  Time  Plot  for 
OFHC  Copper. 
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data.  But  for  the  alight  deviation  beyond  the  second  rebound  of 
the  spall  signal,  the  velocity  vs.  time  plots  compared  well. 

To  demonstrate  the  ability  of  the  failure  model  to  calculate 
the  stress  history  at  the  spall  plane,  we  simulated  a  plate  impact 
test  at  different  velocities.  Various  computer  runs  were  made. 
Figure  49  shows  the  stress  histories  at  V  -  75,  100,  150,  and  300 
m/s.  The  assumed  model  parameters  in  these  runs  were,  s  *»  1  . 0 
(kbar)  oN  =  5.0  (kbar),  B  =  50.0,  f  -  0.01.  It  can  be  seen  from 
these  figures  that  the  reduction  in  stress  due  to  void  growth 
increased  with  increasing  velocities  (or  pressures).  From  a 
threshold  spall  condition  at  V  =  75  m/s,  the  spall  reached  a  well- 
above-  threshold  (excessive  spall  -  the  tensile  stress  relaxed  to 
zero  almost  instantaneously  at  V  -  300  m/s)  condition  at  a  velocity 
of  150  m/ s .  It  is  an  important  feature  of  this  model.  The  simple 
failure  models,  such  as  a  critical  stress  based  model  (i.e., 
failure  occurs  at  o  »  o  ),  will  not  distinguish  the  difference 
between  above  and  below  threshold  spall  conditions.  Models,  in 
which  damage  does  not  degrade  the  strength  (reduction  in  flow 
stress),  are  not  capable  of  modeling  the  dynamic  ductile  failure 
process  realistically. 

The  void  volume  fraction  distribution  in  the  target  plan  is 
shown  in  Figure  50.  The  flyer  plate  is  between  x  -  -3.0  to  0  and 
the  target  plate  is  between  x  =  0  and  x  =  6mm.  The  void  nucleation 
for  low  velocity  impact  is  limited  to  around  1.75  mm  either  side  of 
the  spall  plane  (x  =  3  mm).  The  maximum  f  is  0.0115.  At  higher 
velocity  impacts,  the  nucleation  regions  spread  more  around  the 
spall  plane  and  also  the  void  content,  f  increases  significantly. 
For  instance,  the  f  x  at  V  -  300  m/ s  is  several  times  higher  than 
the  fmax  at  V  -  75  m/s.  For  the  same  model  parameters,  a  change  in 
the  impact  velocity  showed  a  physically  acceptable  changes  in  the 
3tress  history. 

We  showed  that  the  dynamic  failure  model,  based  on  void 
nucleation  and  growth  mechanisms,  could  model  most  of  the  salient 
features  of  spall  signals  obtained  from  plate  impact  tests 
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Figure  49.  Calculated  Stress  History  at  the  Spall  Plane  for  Different 
Impact  Velocities. 
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Figure  50.  Void  Volume  Fraction  (f)  Distribution  in  the  Target  Plate 
at  Different  Impact  Velocities. 


reasonably  well.  We  determined  the  failure  model  parameters  using 
the  VISAR  data  obtained  for  C1008  steel  and  annealed  OFHC  copper. 
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SECTION  V 
CONCLUSIONS 


5.1  SUMMARY 

This  report  described  the  high  strain  rate  results  for 
several  metals.  Data  were  obtained  using  (a)  the  split  Hopkinson 
bar,  (b)  the  long  rod  impact,  and  (c)  the  plate  impact  tests. 

Using  the  induction  coil  heating  technique,  stress-strain  data  at 
high  temperatures  were  obtained  for  C1008  steel,  Nickel  200,  60 6 1  — 
T6  aluminum,  W-2  tungsten,  and  52125  steel.  The  flow  stress 
variation  at  10  percent  strain  with  respect  to  temperature  was 
shown  in  Figure  3*  Several  new  experimental  techniques  that  were 
developed  for  high  strain  rate  measurements  are  described  in  this 
report.  Extended  details  of  the  new  techniques  (such  as  the  man- 
ganin  gauge  technique,  transverse  gauges,  high  pressure 
measurement)  can  be  found  in  the  cited  References  16-30. 

We  successfully  modeled  the  high  strain  rate  data  through  the 
Bodner-Partom  model.  The  B-P  model  parameters  for  several  metals 
were  presented  in  Table  7.  The  spall  failure  was  described  by  a 
state-of-the-art  failure  model.  The  model  was  based  on  microvoid 
nucleation  and  growth  in  materials  due  to  high  triaxial  tensile 
stresses.  The  model  contained  fewer  parameters  (only  four)  than 
most  of  the  other  sophisticated  models  in  which  the  flow  stress  was 
relaxed  due  to  damage.  The  newly  proposed  model  is  continuum 
mechanics  based  and  is  applicable  to  a  general  stress  state  and 
strain  rates.  Preliminary  model  parameters  were  successfully 
determined  for  OFHC  copper  and  C1008  steel. 
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5.2 


RECOMMENDATIONS 


5.2.1  Experiments 

Since  the  exploding  wire  technique  to  expand  circular 
thin  ring/cylinder  geometries  is  perfected  through  several  improve¬ 
ments  in  electronics  and  mechanical  design  and  also  can  be  used 
routinely,  data  for  several  metals  must  be  acquired.  The  radial 
expansion  of  the  ring  is  currently  measured  using  a  streak  camera  - 
flash  lamp  optical  setup.  Preliminary  results  on  1100-0  aluminum 
showed  a  reasonable  streak  record  in  which  the  time  scale  is  well 
resolved.  However,  future  experiments  must  consider  validating  the 
streak  measurements,  through  the  VISAR  measurement  of  the  velocity 
history.  The  long  rod  experiments  must  be  conducted  on  other 
metals.  The  strain  rates  must  be  varied  by  performing  tests  in 
which  the  impact  face  of  the  rod  is  conical  with  the  cone  angle 
varying  between  15°-60°. 

5.2.2  Modeling 

The  long-rod  experiment  must  be  simulated  through  a 
two-dimensional  finite  difference/element  code.  The  validity  of 
the  B-P  model  parameters  must  be  verified.  If  necessary,  the 
initial  yield  stress  measured  from  this  technique  must  be  included 
into  the  B-P  model  parameter  evaluation  technique. 

The  numerical  algorithm  for  the  newly  developed 
failure  model  must  be  improved.  Since  the  void  volume  fraction 
rate  increases  steeply  at  the  finite  difference  grid  points,  it 
introduces  severe  constraints  on  the  time  step.  It  was  also 
noticed,  that  the  results  from  the  current  numerical  scheme  seemed 
to  be  grid  sensitive.  These  numerical  aspects  have  to  be 
investigated.  A  well  defined  failure  model  parameter  evaluation 
scheme  has  to  be  developed.  The  model  must  be  extended  to  other 
metals.  The  generality  of  the  model  parameters  must  be  tested  by 
applying  the  model  to  describe  failure  in  other  specimen  configura¬ 
tions,  such  as  thick  cylinders,  spall  in  a  right  circular  conical 
target  impacted  by  a  flyer  plate  at  the  base,  etc. 
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